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General Asymptotic Expansions of Laplace Integrals 


A. ERDELYI 


1. Introduction 


It is well known [see*, for instance, DoETSscH (1950—1956) vol. 2 Part I, 
WIDDER (1941) Chapter 5, ERDELYI (1956) sec. 22) that, under certain circum- 
stances, an asymptotic expansion of /(¢) as ¢-+0+ induces an asymptotic ex- 
pansion of the Laplace transform 


F(p) a cP" (0) dt (1.1) 


as p->oo, and that a corresponding connection exists between asymptotic ex- 
pansions of /(¢) as ¢—+co on the one hand and of F(f) as # approaches that 
singularity with the largest real part on the other hand. 


The best known expansions in this connection are those involving powers of 
the variables, possibly multiplied by exponential functions. Other expansions 
have also been studied (see the references given above), but on the whole most, 
if not all, asymptotic expansions of Laplace transforms investigated are asymp- 
totic expansions in the sense of PoINCARE: the partial sums of the asymptotic 
expansions are linear combinations (with constant coefficients) of the functions 
of a given asymptotic sequence (scale, gauge). 

Asymptotic expansions of a more general type have been known for some 
time [see, for instance, SCHMIDT (1936) or VAN DER CorPuT (1955—1956)], and 
they have been used occasionally; but they have not been exploited for the 
asymptotic expansion of Laplace integrals. In recent years general asymptotic 
expansions found manifold applications, for instance, in certain singular per- 
turbation problems [see, for instance, KApLUN & LAGERSTROM (1957)| and in 
the study of coefficients of Laurent expansions [WyMAN (1959) and unpublished 
work of the same author]. Their importance seems to be increasing, and it 
appears to be appropriate to apply them to the asymptotic expansion of Laplace 
integrals: this will be done in the present paper. The presentation follows that 
in an earlier note, ERDELYI (1947), on asymptotic expansions in the sense of 
POINCARE, and opportunity will be taken to remove some unnecessary restrictions 
imposed in that note and to correct an error. 

The older results on asymptotic expansions of Laplace integrals have been 
applied successfully in more general methods of asymptotic expansions of func- 
tions defined by integrals, notably in Laplace’s method and the method of 


* References are listed at the end of this paper; they are quoted in the text by 
giving the name of the author followed by the year of publication in parentheses. 
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steepest descents [see, for instance, DoETSCH (1950—1956) vol. 2, chapter 3, 
ErDELyi (1956) chapter II]; and it is not unreasonable to assume that the 
results developed in this paper will permit of similar applications. (See also 
section 9 where an example of this occurs.) 


2. Preliminaries 

Let jf, 2,..<5 9, @,, «9. be functions of a variable z defined on a set R, and 
let z) be a limit point of R (z itself may but need not belong to R). In our case 
z will be either ¢ or #, so that z is a real or complex numerical variable, and it 
is clear what is meant by “limit point’? and “neighborhood” in these cases. In 
the case of ¢, R will be the semi-infinite interval #>0 or t=O, and 2) will be 
either 0 or ov; in the case of ~, R will usually be a region in the complex plane, 
and occasionally a semi-infinite interval on the real axis; in both of these cases 
% will again be 0 or oo. 

We say that {=O (g) in R if there exists a constant (7.e., a number independent 
of z) A so that | f(z)| << A|q(z)| for all z in R; that =O(g) as z 2) if there exists 
a neighborhood U of z) and a constant A so that | f(z)| <A|q(z)| for all z common 
to U and R; and that f=o(y(z)) as zz, if corresponding to each e>0 there 
exists a neighborhood U, of zy so that | f(z)|< «| q(z)| for all z common to U, and R. 

Let {y,} be a finite or infinite sequence of functions @,, defined on R. We 
say that a relation involving the ~,, and possibly other sequences of functions 
or numbers, holds for all n if it holds for all those for which all terms occurring 
in that relation are defined. For instance, if the relation is 29,4;=9,+9,-1, 
and m runs through the sequence of non-negative integers, then “all ”’’ means 
n=1; and if m runs through the integers 0,1,...,.N, then “all »’’ means, 
{<n1<N—1. 

We say that {q,} is an asymptotic sequence if, for all n, ~,.4=0(@,) aS 22. 
When it seems desirable, we say more specifically that {@,} is an asymptotic 
sequence for z—>%); and occasionally we find it useful to indicate the range of n. 

A sequence {y,} is said to dominate {y,} if m ranges over the same set of 
values in the two sequences and if, for all n, g,=O(,) as z>2. These two 
sequences are said to be equivalent if each of them dominates the other. If one 
of two equivalent sequences is an asymptotic sequence, then the other will also 
be an asymptotic sequence. 

The formal series )'f,, is said to be an asymptotic expansion of f with respect 
to the asymptotic sequence {p,} if m ranges over the same values in the series 
>, and the sequence {,} and if for each n, 


f— Lh=0(p,) as &—> &. 


(The sum occurring in this definition is a finite sum so that the left-hand side 
is a well defined function on R.) In this case we write 


(m>) i {p,} aS z—>2. (2.1) 


Sometimes we shall write in greater detail “as %—>2Z in R’’, and occasionally 


we shall omit the indication of the asymptotic sequence or the qualifying phrase 
as 2—>2,y". 
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The expansion (2.1) will be called an asymptotic expansion in the restricted 
sense, or in the sense of Poincaré, if for each n, f, =a, ,,, where the a, are constants. 

The theory of asymptotic expansions in the restricted sense is well known. 
In particular, given f and {@,}, such an expansion, if it exists, is unique, and 
its coefficients can be determined recurrently by Poincaré’s formula 


a, = lim f(z) 4 ba Ap Pr (|/0, (z) 2 
Be AO k<n 

The theory of the more general asymptotic expansions (2.1) resembles in 
many respects the theory of asymptotic expansions in Poincaré’s sense. The 
most conspicuous difference between the two theories is a complete lack of 
uniqueness in the more general case. A function which possesses an asymptotic 
expansion with respect to an asymptotic sequence will clearly possess an infinity 
of asymptotic expansions, in the sense of (2.1), with respect to the same sequence; 
and any asymptotic expansion of a function will be an asymptotic expansion 
with respect to an infinity of asymptotic sequences, for instance with respect 
to all asymptotic sequences which dominate the one originally given. 


Throughout this paper # is a complex variable, and we always write p=o0-+ 10. 
If f(¢) is defined on ¢>0, is integrable (in Lebesgue’s sense) on the interval 
0<t<T for each T>0, and if for a fixed value of p 


is 
F(p) = jim fee j@) di (2.2) 
> 9g 
exists, we say that / possesses a Laplace transform for # and call F(f) the 
Laplace transform of f(t). We shall always denote functions and their Laplace 
transforms by corresponding lower case and capital letters. 

We denote by L the class of all functions f(¢) which possess Laplace trans- 
forms for some # (this # may vary from function to function); by Ly the class 
of the functions which possess Laplace transforms for each =@>0; and by 
Lo the class of those functions in Ly which fail to possess a Laplace transform 
for p=0. 

If f possesses a Laplace transform for fp=0)9+70 9, then it will possess a 
Laplace transform for # whenever o> @,, and F(p) is an analytic function of 
p which is regular in the half-plane > ). It will be useful to recall here some 
other elementary properties of Laplace transforms and to indicate the proofs 
of some of them. 

Suppose that f is in Lo and f(t) 20. Then F(o) exists, for g@>0, not only 
as a limit, as in (2.2), but as a Lebesgue integral (1.1). Moreover, /’(g) is clearly 
a decreasing function of 9. We wish to show that F(g) is unbounded, 1.e., 
F(e) ++ as g->0+ in this case. Indeed, if F(9) <A for e>0, then also 


Fp(o) =f ef) dts A 


for each T>0 o>0. Since F;(g) is a continuous function of g for all real @, 
F,(0)<A. Now F;(0) is an increasing function of T, and since this function 
is bounded (by A), lim F,(0) exists as Tov, so that f possesses a Laplace 
transform for =O and is not in Lo. 

1* 
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Next suppose that f is in L and }()=0. Then F(o) exists as a Lebesgue 
integral for sufficiently large values of @ and isa monotonic decreasing function 
of 9. F(o) +0 as e>+o, but if f()>0 in some neler noe of ¢=0, then 
for each e>0, e°?F (0) >co as g—>-- . Indeed if F(o)sAe ‘ for some APO 
a~0 and all sufficiently large g, then we can use Phragmén’s inversion formula 


[Dortscx (1950—1956) vol. 1, p. 268] 


y oO _A\nh=1 Pan 
f tu) du =1im SF f(n 0) 
a o> 00 A : 
0 


and for 0<¢<a obtain 


; oo 4 ys 
[ fw) dus lim y= e-anoinot 
0 


se) ei n! 


= lim Afexpe *92— 1] =0 
Q@—> oo 


so that f(¢)=0 almost everywhere in 0<t<a. 

For an / in L, the Laplace transform need not converge absolutely for any 
value of #, and the integral in (1.1) may not exist as a Lebesgue integral. Never- 
theless it is possible to transform (2.2), by integration by parts, into a Lebesgue 
integral. Indeed, if F(o) exists for b=fy—O9+70 with o9>0, then 


(pa fe du=o(e*) as too, (2.3) 
and 


FO)=p Let hat ee (2.4) 


where the integral on the right hand side is a Lebesgue integral. 


3. The basic inequalities 


The results on asymptotic expansions can be deduced from certain inequalities 
which are in essence comparison theorems showing that information on the 
relative behavior of two functions at 4=0 or t=oo can be translated into infor- 
mation about the relative behavior of the Laplace transforms of these functions 
at p=oo and at the singularity with the largest real part respectively. The 
two functions will be called g(¢) and h(t) and their respective Laplace transforms, 
G(p) and H(f). 

The first of these inequalities correlates the behavior of an L-function as 
¢—>0 with the behavior of the Laplace transform as f—>oo in a certain region 
in the right half-plane. 


Lemma 1. Suppose that g and h are in L, h(t)=0, and h(t)>0 for 0<t<ht. 
Then 


tim [C@) < fm lel (3.1) 


when t->0+, and p->oo in such a manner that 0—>+0 and for every fixed e>0, 
p=0(0 H(0) e*). 
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Before proving this lemma we note that the last condition is certainly satisfied 
if p is restricted to the sector 


Sa: largp|S3a-—A (3.2) 
with A>0, for in this case |p|/o<cosec4, and H(g) e?—+co as 9p->oo. The 
condition stated in the lemma allows, in some instances, arg p to approach +42. 


Proof. If g and h possess Laplace transforms for 9,>0, then G(p) and H() 
exist, and H(o)>0, whenever o> 0). With a fixed T, 0<T<t#,, we write 
G(p)=1,+1,, where 


T 
Le PO ae, 
0 


A 
Te liml))f ex? oe (0) din 


A>O Tf 


As in (2.4), I, may be converted, by integration by parts, into an absolutely 
convergent integral, 


where 


We set 


and assume that Up<co for some T. (Otherwise the right hand side of (3.1) is 
equal to +00, and there is nothing to prove.) Then | g(¢)|SU;h() for 0<t<T, 
and hence 


ih 
[Zi] feo Uh) dt = Ur H() - 


Moreover, g,(¢) e °’ is bounded when t=T, say |g,(é)|< Be’. If we assume 
0> 20) so that 0 —@)>e, we have 

Smee eu E sore ge hus ice gt 2 le FOL 2.05 
so that 


[Fa] << 221Pt over, 0 > 2Q- 


Thus for any 0< T<#, and o> 20) 


IGP] <p, 4 2BlPl rer 


oe “10 @) 


Now let > in a manner satisfying the last condition of the lemma, and take 
e=+4T. Then the second term on the right hand side approaches zero and 


Since the left hand side here is independent of 7, and the nght hand side tends 
to the right hand side of (3.1) as T—~0-+, this proves the lemma. 
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The second inequality correlates the behavior of f as {-+oo with the behavior 
of F at a singularity with the largest real part. Since F’'(p+c) is the Laplace 
transform of e~“' f(é), we may assume, without restricting the generality of the 
result, that the said singularity occurs at p= 0, and this has been done in formulat- 
ing the following result. 

Lemma 2. Suppose that g is in Ly and h is in Lo, h(t)=0, and h(t)>0 for 
t>t,>0. Then the inequality (3.1) holds as t+ co and p->0 im the half-plane 
o> 0. 


Proof. For any fixed 7 >#, set 


Uy = sup Ol s+ 7], 


and assume that U;<co for some T. (Otherwise there is nothing to prove.) 
The integral defining H(o) for @>0 exists as a Lebesgue integral, and since 


g(t) is dominated by U;h(t) when t=T, the integral defining G(f) for e>0 
also exists as a Lebesgue integral. We again write G(f)—1,+/,, where 

at loo) 

L=fe etd, Gaetan 
0 oe 
For e@>0 we have 
a i 
| SS |e@|dt=Cr, 


say, and 
Lae =) e-'U_ h(i) dt= Up H(0), 
so that 
|G(p)| Cr 
Hig) =e tae 
Now, H(g)->oo as o-+0-+, and hence 
rr Kea) 
lim ae <= U; 


as pO in the half-plane @ >0. Since the left hand side is independent of T, 


and the right hand side approaches the right hand side of (3.1) as T—+co, this 
proves lemma 2. 


The condition h (¢)=0 may be relaxed considerably since in the case of lemma 4 
only small values of t, and in the case of lemma 2 only large values of ¢ matter. 
Indeed, suppose that /, and hy are in L, and h,(t)=h,(t)>0 for 0<t<t#,. Then 


HT, (p) — H,(p) = fe?! [hy (f) — hy (‘)| dt = O(e~*") 
a 
for all sufficiently large @; and for 0<e< t,, we have 
e"* [Hi (@) — Ha(e)] > 0 as e>+o, 


while e** H, (9) and e*® H,(0) are unbounded, so that 


Ai, (@) Ess e°? Hi, ( ) 
AON eet te ee Ht 
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Again, suppose that 4, and A, are in Lo, and hy (t)=h,(t)>0 for t>4,>0. Then 


H(P) — Held) = fe? h(t) — ha ] a 


is an entire function of # which is bounded in some neighborhood of the origin, 
and since H,(0) and H,(0) are unbounded as 9 +0+, we again have H,(0)/H,(0) 1 
as, 0 >-O->: 
4. Asymptotic sequences 
The results of the preceding section may be used to show that under certain 
circumstances the Laplace transforms {@, (p)} of an asymptotic sequence {g,, (¢)} 
form an asymptotic sequence. 


Theorem 1. (i) Suppose that {y,(t)} 1s an asymptotic sequence for t>0+ 
such that for each n, ,, 1s in L, , (t)=0, and ,,(t)>0 for 0<t<t,. Then {®,,(0)} 
1s an asymptotic sequence for p—>oo in such a manner that also 9—>-+-00. (ii) If, 
im addition, there is an unbounded set R in the p-plane such that 9 +-+-co when- 
ever p->oco in R, and if for each n, ®,, (0) =O (®, (p)) as p->0o in R, then {®,, (f)} 
zs also an asymptotic sequence as p—>oo in R, 


Proof. (i) If foo so that g++, we have from lemma 1 


lim —2e+ 


( ) < lim Pn-+1 (2) 
@—> +00 n(Q) ~ t+0+ nit) 


? 


and the right hand side is equal to 0 since {y, ()} is an asymptotic sequence. 
Thus, ®,.,(e)=0(®, (0)) as p->co in this case, and this proves the first half 
of the theorem. 

(ii) Since ,(¢)20, we have in any event |®,(p)|<@,(0), and hence 
®,,(p) =O(®, (0)), as poco in R. Since by assumption also ®, (0) =O (®,,(f)), 
it follows that the two sequences {@, (o)} and {®, (f)} are equivalent as p—oo 
in R. Since the first sequence is asymptotic, the second must also be asymptotic. 


Theorem 2. (i) Suppose that {p,,(t)} is an asymptotic sequence for too such 
that for each n, p, is in Lo, 9, (t)=0, and op, (t)>0 for t>t,. Then {®,,(0)} ts 
an asymptotic sequence for p—>0 in the half-plane @>0. (ii) If, im addition, there 
is a set R in the half-plane 90>0O such that p=0 ts a linut point of R and if for 
each n, ®,, (0) =O(®,,(p)) as p>0 in R, then {®,,(f)} ts also an asymptotic sequence 
for pO im R. 


Proof. From Lemma 2, 


lim PD, +4 (0) ea lim Pn +1 (2) 
e>0+ D,(Q) ~t++00 Qn(t) 


? 


and the proof is quite similar to that of theorem 1. 


5. Asymptotic expansions 
The results of section 3 may also be used to deduce asymptotic expansions 
of Laplace transforms. The technique is well known, and one arrives in this 
manner at results which establish under very general circumstances what 
DoeEtscuH [(1950—1956) vol. 2, p. 43] calls the “ideal case of Abelian asymptotics’, 
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Theorem 3. Let fy, (t)} be an asymptotic sequence, for t>O+, which 1s equi- 
valent to a sequence {p, (t)} satisfying the conditions of theorem ‘1 (i); let R be an 
unbounded set in the p-plane such that 9+ 00, and for each n and each e>0 


2d. ===> 0 (5.4) 
Q P,, (@) ere 


as p—>co in R; and let {X,,(p)} be a sequence equivalent to {D,,(o)} as p>oom R. 
If, under these circumstances, 


tt) ~ dtr) {y, O} (5.2) 
as t>O-+, where f and all the },, ave in L, then 
F(p)~DF,(6) {X,()$ (5.3) 
as p—>ooim R. 
Proof. Fix » and set 
C= fi pe ee (5.4) 
ken 


Then g and / are clearly in L, , (t)=0, and », (f)>0 for 0<t<t?,. Thus, the 
conditions of lemma 1 are satisfied and, moreover, for each e>0, =o (o ®,, (0) e°°) 
by (5.1) as P->coin R. It is a consequence of (5.2) that g(t) =o0(y, (t)) ast>0-+, 
and since {y, (¢)} and {q, ()} are equivalent, also g(t)=o(@,(#)) as t>0+. It 
follows that the right hand side of (3.1) vanishes, and consequently also 
lim G (p)/®,,(0)=0 as p>co in R. But then G(p)=0(@, (e)), and since {®, (0)} 
and {X,()} are equivalent also G(p)=0(X,(f)), as poo in R. Since this 
holds for each n, we have (5.3). Note that {@, (0)! is an asymptotic sequence 
by theorem 1, and hence also {X,,(p)} is an asymptotic sequence *. 

At first the infinite set of conditions (5.1) might seem to restrict R rather 
severely but this is not really the case. As has been pointed out in section 3, 
a sector S,, A>0, is always an admissible R, and in special cases larger regions 
R might be admissible. 


Theorem 4. Let {y,(t)} be an asymptotic sequence, for tcc, which is equi- 
valent to a sequence {,,(t)} satisfying the conditions of theorem 2 (i); let R be a 
set in the half-plane @>0 such that p=0 is a limit point of R; and let {X,,(p)} 
be a sequence equivalent to {®,(e)} as p>0 in R. If, under these circumstances, 
(5.2) holds as t-> oo, and if f and all the f,, are in Ly, then (5.3) holds as p—>0 in R. 


Proof. Fix m and define g and h by (5.4). Clearly g is in Ly and h is in Lo, 
Pn (t) 20, and @, (t)>0 for t>4,. Thus, the conditions of lemma 2 are satisfied; 
and from this point the proof resembles closely the proof of theorem 3. 


6. Sequences of powers 


In the remaining sections we shall give some examples and applications of 
the general results, and will comment on some points arising from the comparison 
of our examples with standard results. 

* The corresponding results in ErpELY1 (1947), theorems 1 and 2, are defective 


in that the sequence denoted there by {y,(p)} need not be an asymptotic sequence. 
The corrected form appears in ErpELYI (1956) p. 31f. 
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The first, and most important, example concerns powers of the variables. 
With complex parameters 2, we set 


Wn (t) as es Pr (t) = fe Pat Re An SSOy 


Clearly, {p, (¢)} and {q, (¢)} are equivalent sequences on the semi-infinite interval 
t>0. The ,, certainly possess Laplace transforms 


co 


Di (pee [ect on dia | (Red) p Re 
0 
for all g@>0, and fail to possess Laplace transforms for f=0, so that all @, are 
in Lo. Moreover, @, (t)>0 for #>0 and all m so that both lemma 1 and lemma 2 
apply here. 

If the ReA, form an increasing (decreasing) sequence, both {y,,(¢)} and 
{p, (t)} are asymptotic sequences for £->0-+ (t->oo), and we can apply theorem 3 
(theorem 4). In view of the special form of ®,, (0), (5.1) is satisfied if and only 
if pe’ *°->0 for each fixed e>0. 

From now on Ry will always stand for any unbounded set in the p-plane 
such that @—>-+oo and for every e>0 pe *®—0, as poo in Ry. We shall 
continue to use the notation Sy for the sector defined in (3.2). 


We also set 
X,(p)=p-™. 


Now, writing =| | ¢’”, e=|/| cos B, we have 


Sak (39 a yy gl 

|Xn (p)| ( ») ore An 3 ( | (cos 9) Re An d 
and see that this quotient is bounded if and only if sec # is bounded, 7.e., if £ 
is restricted to S, with some A>0. Since the above quotient is always bounded 
away from zero, it follows that {Y,(0)} and {X,,(f)} are equivalent in every S4 
with A>0, and in no larger region, as $0 or pow. 


From theorems 3 and 4 we then obtain the following results: 


Theorem 5. For each n, let Red, >0, Red, >ReAd,_,. If 


fO~LtO te} (6.1) 
as t—>0+4, and if f and all the f,, are in L; then 
Cape Sale a Coed (6.2) 
as p->co in some Ry, and 
F(~)~ 2G) {o°*} (6.3) 


as p—>co in some Sy, A>0. 


Theorem 6. For each n let Rei, >0, Reda,<Red,_1. If (0.1) holds as 
t—>-+ co, and if f and all the f,, ave in Ly; then (6.2) holds as p->0 1m the half-plane 
o>0, and (6.3) holds as p—>0 in some S4, A>0. 


Note that the asymptotic expansion F ~ >) F, has been established here in 
two different senses: (6.2) has a larger region of validity, but (6.3) has a better 
asymptotic sequence. This flexibility is lost if one restricts oneself to asymptotic 
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expansions in the sense of Poincaré, for (6.2) can never be such an expansion 


[see (6.5) below]. 
The best known case of these expansions is the particular case in which 


(6.1) assumes the form 
[eyelet Mee, (6.4) 
where the c, are independent of ¢, and (6.2) and (6.3) become 


FQ) Piel UP ie. moe Ime (6.5) 
as the case may be. 
Other special cases of theorem 5 involve inverse factorial series of various 
kinds. With a non-negative integer ”, we have 


—pt eet ndt= n! See 
fe aa: (p+ 1)... (+n) c 
0 
y—bt (pt 4)" dt = n! e) ae 
fe ae p(b—1) ... (b—n) cies 
0 
eee = a Ne AN me (2n)! o>n. 
(pb—n) (p—n+1) ... (b+n) 
0 


Moreover, the three sequences {y, (t)} involved here are equivalent, respectively, 
to {2"},. {2}, ee as ¢->0-+, and their Laplace transforms are correspondingly 
equivalent to {p~""*}, {p-"""}, {p-?"-1 as poco. With constant coefficients, 
c,, we thus obtain the following examples to theorem 5: 

If /(¢) is in L, and if 


KO 2G, (1 a=) {fe Ton te (6.6) 
then 


! n —n— 1 
FU 2 pie i). G-Fo) - oa ae oe 


(inverse factorial series) ; if 


f(t) ~ Dic, (e — 1)” {ml as #304, (6.8) 
then 


N! Cy ae : 
F(t)~> Fe ae {o-""}il'as “pisco in’ R, (6.9) 


(inverse backward factorial series); and if 


f(t) ~ dic, (2sin 24)” {A") ds ie Oe (6.10) 
then 
Ay (2n)! Cy re : 
Me a ee Gia en) Seeds 


(inverse central factorial series). 


. (6.4) is an asymptotic expansion in the sense of Poincaré, and so is (6.5) 
if the second asymptotic sequence is used and accordingly # is restricted to S,, 
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A> 0. The expansions (6.6) to (6.11) may be rewritten as asymptotic expansions 
in the sense of Poincaré by replacing the asymptotic sequences indicated in { } 
by the equivalent sequences formed by the functions occurring in the expansions 
themselves, and again restricting p to S,, A4>0. Thus, these expansions, at 
least in S,,, enjoy uniqueness properties not available in the general case. 


7. Some special Laplace integrals involving logarithms 
For the discussion of asymptotic expansions involving logarithms, we require 
certain results on the Laplace transforms of functions like ¢*~1(-+Llog #)*. Since 
these results are useful in several connections, they will be given in more detail 
than is required for the application in section 8. 


First we discuss the function 


L(o, A, «) J=f(- hoe ue edz (7.1) 


assuming that 0<c<1, A>0, « real; and investigate the behavior of this function 
aS 0—> + oo, 
Differentiating 


aban eh dA) Ge 4 aAS=0, 02 0 


k times with respect to -- A, we first obtain 


foe) k 


f (4log yt? emetdt => (4.1) (2) P(A) e-* (Flog gy’ A>0, e>0 (7.2) 


0 7=0 
where either all upper or all lower signs must be taken. 
Next, we wish to show that 


fo(= log OF Fh are= 2d? S20 ((2.A)* 07 4 Nes, (7.3) 
uniformly in 4 for A>0, as g++ 00. To show this, we remark that |log¢|*< Ad 
for some A>O and all t=c, so that the integral on the left hand side of (7.3) 
is dominated by 


Aufiiver ae 


Furthermore, o¢/(24)<expot/(24) so that the last integral in its turn is domi- 
nated by 


[o.e) 


A(22)" f e-tetde = 2A(22)* e728 ee, 
Q 


a+1 
o*t 
c 


This proves (7.3) since A is clearly independent of A. 
From (7.2) and (7.3) we have for any non-negative integer k, 


L(0,4, 8) => (—1) (*) 1) o-* log e+ O((2aytg~H ee), (7.4) 


v=0 


uniformly in 4, as @—>-+ ©. 
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We wish to extend this result from integer to arbitrary real values of «. 
This extension will be based on the functional equation 


[o.e) 


1 Ya : = bars 
ara L(o,A+2,y) 4z2=L(e,4,y — B), C=O", S06 (7.5) 


0 
which is easily proved by substituting the basic integral (7.1) for L(e, A+-z, y) 
on the left hand side of (7.5) and interchanging the orders of integration. Since 
c<41, the double integral is absolutely convergent and can be evaluated as a 
repeated integral. 


Lemma 3. For A>0, « real, and 9o—>-+ 0, 
Lo, bya) ~S (= 1)" (SPA) o-*log oy" fo * (log g)*—"}. 


Proof. Choose a positive integer k>a«: we shall show that the asymptotic 
expansion stated in lemma 3 holds to k+1 terms. Since k& can be arbitrarily 
large, this proves lemma 3. 


We have «=k—f where k is a positive integer, and B>0. From (7.5), 


é { Z foe) 
L(9, 4,0) =sam-( f+ f )#4*LQA+2 i dz=L,41,, (7.6) 
IN BG, a 
say, where Z is an arbitrary positive number. 
Or HSC< il, 
ae AE Ae diml =i? Z)logt q » =, ae 1, —(z— Z) \loge| J » 
Z Zz 


is a bounded function of ¢ so that 


L, =O f (—loga)tt#¥2-* eet aa) 0—> 0 
0 
and by (7.4), 


L,=O0(9-* 7 (log o)*) go. (7.7) 
We also have 


ae 
1 = 1 1 
ip [#70 (2A+ 24) 9-4 e420) dz = Og 4e-H4) p00, (7.8) 


0 
and hence, substituting (7.4) in L,, 
“ z 


IEP ed -r_ 1 f 4-1) Ay 
DA (7) Cow el” oy f PALO +2) dat 


+ 0(9~’ 7 (log 0)*), i (7-9) 


since the O-term appearing in (7.8) may be included in that appearing in Vets 
Further progress now depends on the asymptotic expansion of the integrals 


Z 
1 en yale Le 
rey f PCL OU+a eras. 
0 
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Apart from the factor @~*, these are Laplace integrals with the parameter log @, 
and their asymptotic expansion can be deduced from theorem 5. From 


gP—-1 Pr) (A+ 2) — ys — {PXASS: A\ O(a actat) 


as z—>0 it follows by theorem 5 that 


Z 


1 pir ni (A ae z) gure. ee dz 


(8) 


ae [rts (A) o~* (log py Ps 2 0(o~* (log oe) 


as @—>-+ cc. We now substitute this expansion in (7.9), reduce the double sum 
to a single sum by means of Vandermonde’s identity 


DH) (= C5?) = (a), 


y+s=n 
and obtain 


19, 4,0) =D (— 1)" (2). @) e~*(log e-" + 0(0~*(log e)?*) (7.10) 


n=0 
as @—>- ov, thus establishing lemma 3. 


Secondly we discuss the function 


loc) 


Mio, A, a) fe dogs) ci" edit, (74) 


Cc 


assuming that « and 4 are real, c>1, o@>0 and 90+. 


In place of (7.3) we have in this case 


k pA—1 ,—ot gy chs 
ii (log #)* #1 e dt =0(— ,) (7.12) 
0 
uniformly for A4>6>0 and @>0. This result follows upon remarking that e~?’<1 
and that given any 6>0, there exists an A>0O, depending on & and 6 and c 


but not on / or g, 


From (7.2), with the lower signs, and (7.12) we have for any non-negative 
integer k, 


M(o, A, k) ->(\r JP (A) 9-7 (= log gk +0 ((A— 0) te?) (7.13) 


uniformly for e>0 and A>06>0. 
The extension of this result to arbitrary real values of « will be based on the 
functional equation 


CO 


1 coh ny os 3 , 
ra M(o,A—2z,y)dz=M(o,4,y—B) B>0, e@>0, (7-44) 


0 
the proof of which is analogous to that of (7.5). 
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Lemma 4. For A>0, « real, and 00+, 
M(o, 7 a) Paps fs a) (A) o*(— log ples fo“ (log ore. 
n=0 “ 


Proof. As before, we set c=k—f where & is a sufficiently large positive 
integer and B>0, and we set 


{ , 2 fee) i 2 = 
Moi he) rg Stl) lek eee 
say, where 0<Z<A. As in the proof of (7.7), 
M, =O f (log s)*#-2-1 e-e* de) 


=0(974*4(— log )*) ee > 0+ 


by (7.13), and with 0< d<A—Z, 
Z 
che PEL O(a = 2) “ae Ot o—->0O+. 
ry | 770 ) SE 
0 
The last two results correspond to (7.7) and (7.8); and from here on the proof 
follows closely that of lemma 3 and need not be given in full. 
For the purposes of the next section, only the first terms of the asymptotic 
expansions given in lemmas 3 and 4, 7.e. the results 


L(@, A, &) 
T'(A) @~* (log @)* 
M(e, 4, «) 
D(A) @~*(—log @)” 


are needed. In both results A> 0 and « is real. 


>14 as o>+oo, (7.15) 


>1 as 9>+0 (7.16) 


8. Expansions involving logarithms 
We first set 
Yn (d) = (— log oye 
Py At) = (log 2)RSEs awe’ Ne 
=0 t>c 


where «, and A,, are complex parameters, and Re A, >0. Under certain assump- 
tions on «, and 4, {y,(t)} and {@, (é)} are equivalent asymptotic sequences for 
t>0-+. Moreover, @, is clearly in L, yy, (t)=0 for ¢>0, and g, (t)>0 for 0<t<c. 
Thus, {@, ({)} satisfies the conditions of theorem 4 (i), and we may apply theorem 3. 
{®, (e)}={L(o, Red, Reaw,)} is then an asymptotic sequence for g++ 00, and 
by (7.15) this asymptotic sequence is equivalent to {(log g)®¢% 9~Re4" which in 
its turn is equivalent to {(log 0)" o~™} for g++ oo. Lastly, one proves exactly 
as In section 6 that the latter sequence is equivalent to {(log p)" p~*"} provided 


that is restricted to the sector S, with some A>0. Thus, from theorem 3 
the following result is obtained. 
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Theorem 7. For each n, let Re 1,,>0 and either Re i, >Re d,,_, and «,, arbit- 
wary or Re A,= Red, <; and Rea,<Rea,.,. If 


th~ > f(t) {(— log)™i™""} as t>04, (8.1) 
and if f and all the f,, are in L; then 


F(p) ~> F,(p) {(log 0)” o~%} as p—+oco insome R,, (8.2) 
and 


F(p) ~ > F, (p) {(log f)*™p-™} as p—-ocoin S,, A>o. (8.3) 
Similarly, we may start with . 


Yn (2) 
Pn (t) 


I 


(log £)% ¢4»—4 


| 


I 


(log [Reem fRe An—1 iS oS 4 
O ORIGGC, 


where Re /,>0, and note that under certain assumptions on «, and /,,, {y, (é)} 
and {y, (¢)} are equivalent asymptotic sequences for t+ ce, and {g, (é)} satisfies 
the conditions of theorem 2. Thus, {®, (g)} is an asymptotic sequence for p—>0 
in the half-plane @>0, and (7.16) in connection with an argument similar to 
that preceding the statement of theorem 5 may be used to replace {®, (@)} by 
{(— log 0) o~*"} in Ry and by {(— log p)™p~*} in Sy. Finally one obtains 

Theorem 8. For each n let Red, >0 and etther Rei, <Red,_, and u, 
arbitrary or Reij,=—Red,_, and Rew, < Rew,_1. If 


Ti 


fO ~ Dita. -tlogiris 7}. as. i> 4.00, (3.4) 

and if f and all the f, are in Ly; then 
p) oy Dirt) ie ls 108-0) 07 72} (3.5) 

as p->0 in the half-plane o> 0, and 
Ep) igh enigake pepe} (8.6) 


asp—>Oim Sy, A>0. 

These results appear to be new even in the restricted case when (8.1) and 
(8.4) are asymptotic expansions in the sense of Poincaré. If the «, are not integers, 
an interesting situation arises in this case. If we have, say, 


~ > a,(— logit" as t>0+ 


in the sense of Poincaré, with constant a,, we obtain, with an extension of (7.1) 
to complex «, A, and #, 


b) ~ >a, LP Ans tn) —« {(log e)*" 0} 


as p—>co in Ry. Here the terms of the expansion depend on the choice of ¢ 
used in the construction of L, and in point of fact a different c, could be used 
for each n. Thus, in the case of general «,, the asymptotic expansion of F() 
is not unique, and it is impossible to single out any particular expansion as 
being simpler or better than the others. Moreover, one has the choice of regarding 
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the asymptotic expansion of F(f) either as an asymptotic expansion in the 
sense of Poincaré with respect to the comparatively involved and non-unique 
asymptotic sequence {L(p, A, %,)} for p > in S,, or else asa general asymptotic 
expansion with respect to the comparatively simple asymptotic sequence 
{(— log 0) 9} and valid in the larger region Ry. There is little doubt about 
the second point of view being preferable in most problems. 

The situation simplifies considerably if all the «,, are integers so that (— log 1) 
possesses a unique determination for all ¢>0, and c=oo may be taken in (7.1) 
for all n. If, moreover, the «, are non-negative integers, the expansion (7.2) 
may be used, and the asymptotic expansion of F(p) in S, becomes either an 
asymptotic expansion in the sense of Poincaré with respect to the asymptotic 
sequence 


an 
fS' (—1' (4) P(A,) p-* (log p)*-*} 
r=0 
or else, upon rearrangement of the terms, an asymptotic expansion in Poincaré’s 
sense with respect to the sequence {p “"(— log p)*"~’} in which m runs through 
the same values as in (8.1), for each 2, y runs through 0, 1, ..., «,, and the terms 
of the sequence have been rearranged, if necessary, in decreasing orders of 
magnitude as poo, and eliminating duplications. 

In the literature known to the present writer, only the case of non-negative 
integer «, has been considered, and precise results are known mostly when 
%,=0 or 1 [see DoETscH (1950—1956) vol. 2,-p. 45—50]. 


9. Some Laplace integrals involving exponential functions 
For the discussion of asymptotic expansions involving exponential functions 
we require some information on the behavior of the Laplace transforms of 
functions like ¢’~*exp(«#t’) as 0 or pce. The results which we require 
are not essentially new but there does not seem to be a convenient reference 
giving all of them; and we propose to prove them here, especially since their 
proof illustrates the application of the results of section 6. In point of fact, 
we use the results on Laplace integrals to provide a justification of Laplace’s 
method [DoETscH (1950—1956) vol. 2, p. 83£; ERDELyI (1956) sec. 2.4. or 

the sake of simplicity we restrict ourselves to real parameters and set 


h(t, «, B, 4) =t*—" exp (at) (9.1) 


H(o, «, B, A) car t}—* exp (a2? — 0 #) dt, (9.2) 


the latter provided the integral exists. 


We shall use # as a comparison function, and it is necessary to satisfy the 
conditions of lemmas 1 and 2. In the case of lemma 2, too. In order to make 
h differ in its asymptotic behavior from ¢*~!, we must assume B>O and a=20¢ 
integrability at #=0 demands A> 0, and the condition that h possess no Laplace 
transform for =0 clearly demands «+0. Since «>0, h will possess a Laplace 
transform for p=o>0 only if B<1. In the case of lemma1, t-+0+, and the 
asymptotic behavior of A will differ from that of t4-1 if B<0 and #=t0. Inte- 
grability at #=0 then demands «=< 0, and no condition need be imposed on A. 
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Accordingly, we shall investigate the asymptotic behavior of H(g, «, B, A) 
under two sets of conditions: (i) «>0, 0<B<1, A>0, and 90+; and (ii) 
“<0, B<O, 4 arbitrary (real), and @—oo. In either of these cases the integral 
in (9.2) clearly converges for all o> 0. 


In both cases we introduce the notations 


TAO (9.3) 


note that y and z are positive, and z->+ oc in both cases. Let us introduce a 
new variable of integration x in (9.2) by the substitution t= (« 8 9)’ x, obtaining 
00 
H(o,0,8,2) = (2) fe rerse PH de, (0.4) 
0 
We wish to investigate the behavior of this integral as z->oo. 


Now, the integral in (9.4) is not a Laplace integral, and so theorem 5 does 
not apply directly. Nevertheless, it is possible to transform this integral into 
a Laplace integral by the substitution 


yaa Pra —14 64 (0:5 
and then apply theorem 5. 

The function of x on the right hand side of (9.5) is strictly decreasing when 
x<1 and strictly increasing when *>1, and it becomes necessary to break the 
integral in (9.4) in two parts, and transform each part separately. 

The substitution (9.5) maps the interval 1< *< o onto the interval 0< y< 0, 
and we have 


[oe 


(Lee ee BaP) dy, — 9 (L-B Shly)e**dy 


where 


(9.5) also maps the interval 0< x<1 onto the interval B+ —1>y>0 if 0<f<1, 
and o>y>0 if B<0, so that 


1 foe) 
fat Be PP dy = — OP f fy(y) eM dy, 
0 0 


where 


fly)=e 5. O<*<1, 


andy s(y)\—= Oi for 7p 4— 4 incase (i): 
Setting /(v) =f, (v) —f2(y), we thus obtain 
H(0, 0, B, 2) = (SPY "2-0-0 f f(y) en? dy, (9.6) 


a representation of a multiple of H by a Laplace integral in which z+ oo in 
both cases (i) and (ii). It is clear from (9.4) that f is in L, and it follows from 
Arch, Rational Mech. Anal., Vol. 7 2 
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theorem 5 that an asymptotic expansion of f in series of ascending powers of y 
will lead to an asymptotic expansion of the Laplace integral in descending powers 
of z. Only the leading term of this asymptotic expansion will be computed here 


explicitly. 
Lemma 5. I} either (i) «>0, O<B<1, A>0 and 90+ or (ii) «<0, B<O, 
A arbitrary, and 9-—>-+ 0, then 


H(o, t, B, A) ~ (2B ewmry (ee Cy go rye (9.7) 


where y, z are defined in (9.3), the c, are independent of z, and depend only on B 
and A, Co= 24, and (9.7) is an asymptotic expansion in Poincaré’s sense. In 
particular 


(9, a, B, A) oh 4 9.8 
Va y(a py"? 7 exp— (1B) 2 a 


under either of the conditions (1) or (ii). 


Proof. The right hand side of (9.5) is regular at x=1, the point corresponding 
to y=0. Using (w) as a generic symbol for a power series in w such that 
8 (0) =0, we have by explicit computation 


y=4F (x —4)® 1+ Be — 1) |x —1| <1. 
By inversion of power series we have, in some neighborhood of y=0 
+ (* —1) =)2yy 1+ 8(Vy)], 


where |/y denotes the positive square root, the upper sign is to be used for x>14 
(z.e., in f,) and the lower sign, for O0<x<1 (i.e., in f,). Hence 


ax 


| 4+ ROM) 


(2 


and since 


v4 = 414 B(x — 1) =14+ BV), 


we finally obtain 


say, where Co=V2y, and the power series converges, and represents /(y), in some 
punctured neighborhood of y=0. 

(‘Under these circumstances })c, y"~"? represents f(y) asymptotically, in 
Poincaré’s sense, as y->0+, and it follows from theorem 5, or (6.5), that 


[tolerrtdy~ ye, P(E se as 2 > 400, 
0 n=0 


again in Poincaré’s sense. Substitution in (9.6) leads to (9.7), and (9.8) is a 
simple consequence of (9.7). 

It may be mentioned that certain well-known asymptotic expansions of 
modified Bessel functions and of parabolic cylinder functions are special instances 
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(for B=—1 and B=} respectively) of (9.7); and that the validity of (9.7) can 
be extended easily to complex values of the parameters. In particular, (9.7) 
is valid for complex A with Re A> 0 in case (i), and no restriction on 4 in case (ii). 


10. Expansions involving exponential functions 
By taking appropriate sequences of «, 8, A in (9.1) we can now construct 
asymptotic sequences, in case (i) for {+c and in case (ii) for +0. It follows 
from lemma 5 that these sequences satisfy the conditions of theorems 1 and 2, 
since h is clearly in Lo in case (i) of lemma 5, and in L in case (ii), and h(t) >0 
for #>0. Thus, the Laplace transforms also form an asymptotic sequence. By 
(9.8), the latter sequence is equivalent to the corresponding sequence of functions 


X (0) = oP Ay e(B7—l)z 


as @->0-+ in case (i) and g++ co in case (ii). An examination of the exponential 
part shows that for complex #0 in case (i), and —>co in case (1i), X(0)/X(p) 
cannot remain bounded unless* arg #0 so that there is no useful region R 
in the complex plane in which the sequence of the X(p) is equivalent to the 
sequence of the X(o). 

For the sake of simplicity, we restrict « and f to real values, allowing / to 
be complex. 

In case (i), we set 


yw, (t) = t—1 exp («,, 2) , 
Yn (1 = Ba) s (10.1) 
Xe (0) =0— ™ exo lB, — 1) (On PaO) |» 
and verify that under the conditions given in theorem 9, {y, (¢)} is an asymptotic 


sequence as t->oo, {X,,(0)} is an asymptotic sequence as g->0-++, and the con- 
ditions of theorem 4 are met. Thus we obtain 


Theorem 9. For each n let (i) «,>0, 0<B,<1, Red,>0; and either (ii) 
Ba Pu—19 or (iii) Bn=Bn-1 and Ky by-1> or else (iv) Br=Bn—1 Ay, —1> and 
eA <= REAL. Lf 


1O~XLhE) als (10.2) 
as t->-+ 00, where y,, is given in (10.1), and tf f and all the f,, ave in Lo; then 
F(p)~L Ef) {Xn (0)} (10.3) 


as p->0 in the half-plane 0>0, X,, being defined in (10.1). 
In case (ii) we change the notation slightly, setting 
Pn (t) = te? exp (— a, £-P») 
Ye Pale (10.4) 
X,, (g) = @7 tn exp[— (Bai + 1) (Bn orm)” - 


x In a special case this examination was carried out by DorErscu (1950—1956), 
vol. 2, p. 94. 
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Under the conditions set out in theorem 10, {y,,} is an asymptotic sequence as 
t+0-+, {X,,} is an asymptotic sequence as goo, and the conditions of theorem 3 


are met. We then have 


Theorem 10. For each n let (i) «,>0, B,>0; and either (ii) B, >By—1, or (i) 
Bu=Bn—-1 and «>, 1, or else (iv) B,=By—1, %=%-1, and Red, >Red,_4. 
If, under these circumstances, and with wp, given in (10.4), (10.2) holds as t>0-+, 
and if f and all the },, are in L; then (10.3) holds, with X,, as im (10.4), as poe 
in some Ry. 

This is a straightforward deduction from theorem 3. Since for each m and 
each e>0, X,,(@) is bounded below by some multiple of e~?°?, condition (5.1) 
is Clearly satisfied in Ro. 


The writing of this paper was sponsored by the United States Army under Contract 
No. DA-11-022-ORD-2059, Mathematics Research Center, United States Army, 
Madison, Wisconsin. 
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Uber den Resonanzbegriff bei Systemen 
von n linearen gewobnlichen Differentialgleichungen 


erster Ordnung 


RAHMI IBRAHIM IBRAHIM ABDEL KARIM 


Communicated by L. COLLATZ 


§ 1. Problemstellung, Hauptergebnisse 
R. IGLiscu untersuchte den Resonanzbegriff bei linearen gewéhnlichen Diffe- 
rentialgleichungen zweiter Ordnung (Arch. Rational Mech. Anal. 3, 179—186 
(1959)). Diese Betrachtungen sollen ausgedehnt werden auf Systeme von n ge- 
wohnlichen linearen Differentialgleichungen erster Ordnung, die in Matrixform 
mit x(¢) als gesuchtem Vektor geschrieben werden sollen: 


(1) =U re +O; 


dabei soll die quadratische (z.B. reelle) Matrix (¢) etwa stetig in ¢ und 
periodisch mit der Periode P sein, d.h. alle n? Elemente a;, (¢) sind stetige mit P 
periodische Funktionen in ¢; der (z.B. auch reelle) Vektor j(¢) sei etwa gleich- 
falls stetig und mit P periodisch: 


(2) We+P)=—=AM); feE+tP)=TO. 
Das zu (1) gehdrende homogene System ist 
a) 9 
(3) a 
und das ,,adjungierte“ 
d 
(4) Se ae 


wobei durch das hochgestellte T der Ubergang zur transponierten Matrix ge- 
kennzeichnet sei. 

Definition 1. Bei dem inhomogenen Differentialgleichungssystem (1) liegt der 
Resonanzfall vor, wenn das adjungierte System (4) mindestens einen mit P 
periodischen Lésungsvektor 3 (¢) besitzt, fiir den 


P 
(5) Ss Of Odt=C +0 
0 
ist. 
In §4 wird bewiesen: 


Satz 1. Im Resonanzfall nimmt jeder Losungsvektor x(t) von (1) mit unbe- 
schrankt wachsendem t beliebig groBe Betrige an. 


Ze Raumi IBRAHIM IBRAHIM ABDEL KaRIM: 


Definition 2. Bei (1) liegt der Hauptfall vor, wenn (4) keinen mit P periodi- 
schen Lésungsvektor 3 (t) besitzt. 

In §5 ergibt sich: 

Satz 2. Im Hauptfall gibt es fiir alle t beschrdnkt bleibende Losungen r (¢) 
von (1), 2B. die eindeutig existierende mit P periodische Lésung. 

Definition 3. Bei (1) liegt der Ausnahmefall vor, wenn zwar (4) mit P periodi- 
sche Lésungsvektoren 3(¢) besitzt, jedoch fiir alle diese mit P periodischen 
Lésungen gilt: 

(6) fa OO at=0. 


0 

In §6 wird gezeigt: 

Satz 3. Auch im Ausnahmefall gibt es fiir alle Werte von t beschrankt bletbende 
Lésungen x(t) von (1), z.B. die mit P periodischen, die aber jetzt nicht mehr ein- 
deutig bestimmt sind. 

Die §§ 2 und 3 sind Hilfsbetrachtungen iiber das inhomogene System (1), 
vor allem aber iiber den Zusammenhang zwischen dem homogenen System (3) 
und dem adjungierten System (4) vorbehalten. Dabei wird die Periodizitats- 
voraussetzung (2) erst in §3 eingefiihrt. 

Viele Ergebnisse sind bekannt, werden aber hier auf eine Weise abgeleitet, 
die keine besonderen Vorkenntnisse verlangt. 


§ 2. Hilfsbetrachtungen tiber Systeme 
linearer Differentialgleichungen erster Ordnung 


In diesem Paragraphen wird die Periodizitatsvoraussetzung (2) nicht benutzt. 


0, (), --., 0, (4) sei an der Stelle ¢ ein linear unabhadngiges Lésungssystem 
(Fundamentalsystem) von (3), das wir zu der Lésungsmatrix 


(7) Y (4) = (1, ---, Dn O)) 

zusammenfassen wollen. Dann ist 

(8) Det 2) (to) + 0. 
Satz 4. Aus (8) folgt 

(9) Y@) = Det (¢%) +0 

fiir alle Werte von t. 


Beweis. Bezeichnet man die Adjunkten (Unterdeterminanten mit richtigem 
Vorzeichen (—1)'*”) zum Element y;,() in (7) oder (9) mit Y;, (¢ (t) und setzt 


(10) su={t fiir pak. 
4 Shirt es hs 


so rechnet sich unter Beachtung von (3), wenn die Differentiation nach ¢ durch 
einen Strich gekennzeichnet wird, 


dy 
dt. = 2%, = 24:8 Yaw Yi iv = 240 dee Y = Yas 


v, 1, k 
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und somit 


\ 


Y(t) = Y (4) exp/ i (411 + do2 + +++ + ayy) dt). 


0 


Hieraus folgt unmittelbar Satz 4. 
Satz 5. Neben (7) ist auch 
(11) Y*() = YA) -C 
bet beliebiger konstanter Matrix © mit von Null verschiedener Determinante C ein 
Fundamentalsystem von Lésungen von (3). 
Beweis. (7) und (3) kann man mit 
(12) DY! = (01, Yar ---> De) 


zusammenfassen zu 


Ch ney hee 
(13) SG = UY. 


Damit rechnet sich 


(YC) =YCH=AYC=A(YO). 
Jedes Y)€ ist also Lésungssystem von (3). 
DaB auBerdem det (2)C) =£0 ist, folgt aus 
det (©) = Y-C =0, 
Satz 6. Die Vektoren 


Y,; () 
Y,, (¢ : 

(14) Nea, ai) ; U— AR et; 
Y, ;(0) 


die man in die Matrix 
(15) 8) = (pay Ya) 
zusammenfassen kann, bilden ein Fundamentalsystem von Lésungen von (4). 


Beweis. Zunachst sieht man sofort an (15), daB 
(16) g°- Y=, 


d.h. die Einheitsmatrix ist. Damit kann man unter Einfiihrung der reziproken 
Matrix Y)+ auch schreiben 


(17) B7=9, Br =(97)’. 
Aus (16) oder (17) folgt unmittelbar: 
Zo= eu, 0), 


DaB 8 Lésungsmatrix von (4) ist, folgt so: 
GemaB (16) ist 
8° N+ BY’ =0, 
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also mit (13) _arg=3ray, — ar — stag yr=ary, 


mithin 
(18) Boz = 18, 
was mit (4) tibereinstimmt. 
Aus den Sitzen (5) und (6) folgt nun sofort 
Satz 7. Allgemeines Fundamentalsystem von Losungen von (4) tst die Matrix 
(19) 3*=3- C=C 


bei beliebiger konstanter Matrix © mit von Null verschiedener Determinante C. 


Bezeichnet man 


(20) L(y) = — XW 
und 
(21) L*(3) =— BU (3, 


so gilt fiir zwei beliebige Vektoren y (¢) und 3 () 


d 
(22) 3? L(y) — (£*@))?-9= 4G" -») 
und 

d 
(23) (L(y))"-4— 9" L*@) = - 4). 
Daher gilt die (erweiterte) Lagrangesche Identitat 

d d 

(24) a7 6 19) Os aL. 


fiir Lésungsvektoren t) und 3 von (3), bzw. (4). 
Analog gilt fiir Losungsvektoren von (1) und (4) 
d 


(25) Oo \ad eee oe elena 


§ 3. Uber Systeme linearer Differentialgleichungen 
mit periodischen Koeffizienten 


In diesem Paragraphen ist die Periodizitaétsvoraussetzung (2) wesentlich. 


Satz 8. Die homogenen Systeme (3) und (4) haben die gleiche Anzahl 0 linear 
unabhangiger mit P periodischer Lésungen (0So<n). 


Beweis'. @ sei die Anzahl der mit P periodischen linear unabhingigen 
Lésungen von (3). Es braucht nur gezeigt zu werden, daB (4) genau o linear 
unabhangige Lésungen besitzt, die mit P periodisch sind. Denn der Schlu8 
laBt sich umkehren, da (3) das adjungierte System zu (4) ist. 


Fall 1. Es sei g=n. Die ganze Matrix 9)(¢) ist hier mit P periodisch, daher 
auch die Matrix 3 (¢) nach (15) 


LS VionwRealeniscrs 
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Fall 2. Es sei 0<o<n. Seien die y,(é) mit y=1, 2,...,@ mit P periodisch. 
Dann gilt fiir jede beliebige Stelle ¢, 


(26) yp. (@) et? =o,(44+ P) —», (4) +0 fir w=e+41,...,n. 


Stande namlich in (26) das Gleichheitszeichen, so wiirde aus (3) unter Beachtung 
von (2) die Periodizitat dieses 4,,(f) mit der Periode P folgen. 

Im folgenden sei ¢, eine beliebig gewahlte, dann aber festgehaltene Stelle. 
Integration der ersten Gleichung (24) zwischen ¢, und ¢,+ P liefert fiir jedes 
NLD ors, 


(27) pe) Bk Geli OM Gibray =113;2.6k7, 0% 


Das sind m lineare homogene Gleichungen, die mindestens die @ linear unab- 
hangigen Lésungen , (¢,) fir y=1, 2,..., 9 besitzen. Es soll jetzt gezeigt werden, 
daf} kein weiterer davon linear unabhangiger Loésungsvektor u, existieren kann. 
Sei im Gegenteil angenommen, daf fiir ein solches 1, gleichfalls 


(28) se O|n** uy, =0 
gilt. Dann bestimme man einen Vektor u(¢) aus (3) 
(29) 2 Ui mit al (ei Uy « 


Fiir diesen Vektor u(t) folgt aus der ersten Gleichung (24) durch Integration 
zwischen ¢, und ¢,+ P fiir k=1, 2,...,” 


an (4 +P) u(t) + P) — 32 (4) u (4) =0 
und nach Subtraktion von (28) 
3z (t, + P) [u (4, + bb u(t) | =0. 


Das ist ein lineares homogenes Gleichungssystem mit von Null verschiedener 
Det 1 te; es muB also 
eterminante; es mu CR rae 


sein, und daraus folgt vermége (29) und (2) die Periodizitat von u(t) mit P. 
Wegen (26) ist also 1;=u(t,) von vj, (4), ..-, 9. (4) linear abhangig. 

Nun wissen wir also, daB das Gleichungssystem (27) genau o linear unab- 
hangige Lésungen besitzt. Daher hat die Matrix 


(CV ARS (Gs (2) RSS, Oe oe. 2) 


den Rang ~—o. Durch entsprechende Numerierung ]aBt sich erreichen, daB 
gerade 3,(¢)|2 1°”, ..., 3, —o(é)|2*” linear unabhangig sind. Die tibrigen dieser Vektoren 
kann man durch diese linear ausdriicken: 


gaa) 


n—0 
(30) Be ee Gases tL Bm 0A, M 


Nun betrachte man die Vektoren 


jm (¢) fr 47 = 1,2, 6..;% — 0 


1 3m (t) = n—o 
ee 3m (t) —V>" Cn, obo (2) fir. m=n—o-+1,...,n. 
o=—I 
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Da 
Det (3#(é), .--.4%(@)) = Det (1 0), ++-. 8x) 0 


ist, sind die » Vektoren (31) linear unabhangige Losungen von (4) fiir alle ¢ 
(vgl. dazu auch Satz 4). Die 3,°(¢) mit y=1, 2,...,.%—@ sind wegen der linearen 
Unabhiangigkeit der Vektoren 


einschlieBlich ihrer Linearkombinationen nicht mit P periodisch. Dagegen be- 
sitzen die 3*() mit w=n—o+1,...,m gerade die Periode P, wie aus (30), (4) 
und (2) abzulesen ist. Damit ist unser Beweis im Fall 2 erbracht. 


Fail 3. 9=0. Hatte Gleichung (4) eine mit P periodische Lésung, so miiBte 
nach dem Friiheren auch Gleichung (3) eine solche besitzen. 


§ 4. Der Resonanzfall 
Beweis von Satz 1. Es sei also 3(¢) ein mit P periodischer Lésungsvektor 
von (4), fiir den (5) gilt. 
Wir machen im Gegensatz zur Behauptung die 
(32) Annahme: |x(t)| SE 


fiir alle ¢ und eine beliebig herausgegriffene Losung r(¢t) von (1). Integration 
der ersten Gleichung (25) zwischen ¢ und ¢+mP mit willkiirlichem positiv 
ganzem m liefert unter Beachtung von (5) 


(33) a) [z+ mP)—rQ]=m-|Cl. 
Infolge der Periodizitat von 3’ (¢) gilt fiir alle ¢ eine Beschrankung der Form 
(34) Parry iesa ey 
Damit ergibt sich aus (33) eine Abschatzung 
Ds2B= m|Olyq ee BIEL, 
Das liefert bei geniigend groBem m einen ae ante zu (32). 
Man kann das Ergebnis auch so aussprechen: 


Satz 9. Im Resonanzfall gibt es fiir jeden Losungsvektor x(t) von (1) in jedem 
Intervall 


(35) tp t<St,tmP 
mindestens eine Stelle t*, fiir die 

6 *) > m|Cl 
(36) r)275 


ausfalit mit C aus (5) und D aus (34). 


§ 5. Der Hauptfall 


In diesem Paragraphen wird vorausgesetzt, daB Gleichung (4) keine mit P 
periodische Lésung 3(¢) besitzt. Zundchst beweisen wir 
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Satz 10. Sind 3, (t), ..., 3, (¢) linear unabhangige Lésungen von (4), so ist 
at (a)|,*? 
32 (t)|r” 


(37) D (t) = Det +0 fiir alle t. 


n (a) 


Beweis. Ware fiir einen speziellen Wert t=¢, im Gegensatz zur Behauptung 
D (tz) =0, so kénnte man m Zahlen a, a, ..., %, mit af+---+«2>0 so bestimmen, 
daB 


Di ty iy (2+ P)= 2% o (ty) 


ist. Setzt man 


so ware 


Hieraus wiirde aber nach (4) und (2) 
30 (¢ + P) = 30 (4) 


fiir alle ¢ folgen im Widerspruch zur Voraussetzung. 


Beweis von Satz 2. Ware r(t) eine mit P periodische Lésung von (1), so 
wiirde Integration der ersten Gleichung (25) zwischen einer festen Stelle ¢, und 
t,-+ P liefern: 


ts+P 


(38) alent we = Sir Oty fit 1p. A 


Da die Koeffizientendeterminante dieses Gleichungssystems nach Satz 10 von 
Null verschieden ist, la8t sich hieraus y(¢,) eindeutig bestimmen. Wir denken 
uns nun (1) mit diesen Anfangswerten an der Stelle ¢; gelést und erhalten so 
einen Lésungsvektor x(t). Wir brauchen nur noch zu beweisen: Dieses yx (¢) ist 
mit P periodisch. Dazu ist in Hinblick auf (2) nur ndtig, festzustellen, dab 
von selbst 


(39) t(ts + P) =2 (4) 
ist. Integration der ersten Gleichung (25) zwischen ?; und /,-+P liefert fiir 


DeSSsPs preg 
t+P 


ay (tg +P) - 4 (t+ P) — ap (ts) ° = Sarl 
Zieht man hiervon (38) ab, so erscheint das lineare Gleichungssystem 
(40) gy (ts ++ P) [x (2 + P) — 2 (4)] =0 


mit von Null verschiedener Determinante analog zu (8). Also folgt das Ver- 
schwinden der eckigen Klammer und damit (39). 


28 Karim: Resonanz bei linearen Differentialsystemen 


§ 6. Der Ausnahmefall 


Hat (4) genau 9 linear unabhangige mit P periodische Lésungen 3 (2), ..., 3 (4) 


fiir die 
t+P 


(41) fot dr= 0; 7 =A) 2, 250 


gilt, so hat die Determinante (37) genau den Rang »—g. Von dem Gleichungs- 
system (38) sind die ersten @ Gleichungen (fiir »=1, 2, ..., @) von selbst erfillt, 
da nur Nullen darin vorkommen. Aus den restlichen Gleichungen (38) mit 


y=o+1,...,” kann man g linear unabhangige Vektoren 

(42) UL (ts) » U3 (ts), «++» Le (ts) 

bestimmen, dazu durch Lésung von (1) mit den entsprechenden Anfangswerten (42) 
(43) rr (t), 22 (1), «++, Ee). 

DaB fiir jedes dieser x(t) (u=1, 2, ..., @) von selbst 

(44) tn(ts + P) = ri (6) 

wird, folgt aus dem analog (40) abzuleitenden Gleichungssystem 

(45) ts (to +P) [tilts ++ P) — ur(ts)] =0 


mit von Null verschiedener Determinante. Alle diese x*(¢) haben also die Periode P. 
Wir koénnen somit feststellen: 


Satz 11. Hat (4) genau o linear unabhingige mit P periodische Lésungen, 
fiir deren jede (6) gilt, so besitzt Gleichung (1) eine o-parametrige mit P periodische 
Lésungsschar. 

Dies Ergebnis steht in Einklang mit der trivialen Tatsache, da8 man alle 
mit P periodischen Lésungsvektoren von (1) erhalt, indem man zu einem von 
ihnen sémtliche mit P periodischen Lésungsvektoren von (3) hinzuaddiert. 
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Uber den Resonanz fall bei S  ystemen 
von n nichtlinearen gewohnlichen Differentialgleichungen 
erster Ordnung 


RAHMI IBRAHIM IBRAHIM ABDEL KARIM 


Vorgelegt von L. COLLATZ 


§ 1. Problemstellung, Hauptergebnisse 
In Erweiterung von Untersuchungen von R. Icriscu [1] wird das nichtlineare 
Differentialgleichungssystem 


(1) Eee, (2 Saray re) (Cae 


betrachtet, das sich unter Einfiihrung der Vektoren 


Uy h, (¢) 27 (Bt ys b) 
(2) leo 2 = Tis?) und g(u, 7) = Balm, 22” y 
Uy h, (0) Sy (Uys. 24 My 2) 
in der Form 
(3) Ze = q(u, t) +h) 


schreiben laBt. Es gelte die Periodizitatsvoraussetzung 
(4) g(ut+P)=g(u,4), b¢+P)=5(). 


Fiir die Funktionen g; sei hinsichtlich der Veranderlichen uw; eine Taylor-Ent- 
wicklung bis zu den Gliedern zweiter Ordnung moglich, {(¢) se1 etwa stetig. 


Es sei eine mit P periodische Lésung 1 (¢) von (3) bekannt: 


(5) M0 = g (up, t) +5 (),  tolt-+ P) = 1) (0). 


Untersucht werden dann bei geringer Abanderung von f(é) die zu u(t) benach- 
barten Lésungsvektoren 1 (¢) des Gleichungssystems 


(6) = 9 (u,) +50 +640 


bei gentigend kleinem ||; dabei soll wieder gelten: 


(7) 1e+P) =F. 
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Mit 

(8) u(t) = p(t) + 2 (4) 

geht (6) nach Subtraktion von (5) iiber in 

(9) 7 — 4 (g(t) +4 0.4) —9(tol,4) +81. 


Untersucht werden hiervon die Lésungen mit kleinem |r(¢)| bei kleinem |]. 


Wie in §2, wo eine Umformung von (9) vorgenommen werden soll, gezeigt 
werden wird, spielt die Rolle des homogenen linearen Gleichungssystems* das 


System 
(10) 1) =X) n() 


mit der Matrix'(?, R=4, 2, sy.%) 


0» 0g (Uo, ¢ Og (Uo, t 
(11) 2 () = (aix(¢)) = (P8Re-D) (ee t) : G (Uo, 2) yoo, 20M), 


OUR Ou, OU: OUn 


Offenbar ist 
(12) W(é¢+ P) =A(A). 


Das zu (10) adjungierte homogene System ist dann, wenn durch das hoch- 
gestellte T der Ubergang zur transponierten Matrix gekennzeichnet wird: 


B= a5. 


at 


(13) 


Satz 1 (Resonanzfall). Besitzt (13) einen mit P periodischen Lésungsvektor 
3(t), fiir den 


t+P 


(14) Se QF dt=C +0 


ausfallt, so nimmt jede Lésung x(t) von (9) — unabhdngig von den Anfangswerten — 
mit wachsendem t Werte an, deren Absolutbetrage mindestens von der GréBenordnung 
/|B| sind. — Der Beweis wird in §3 erbracht. 


In §4 ergibt sich dann 


Satz 2 (Hauptfall). Hat (13) keinen mit P periodischen Lésungsvektor 3 (t), 
so besttat (9) zw jedem hinreichend kleinen B Lésungsvektoren x (t), die fiir alle Werte 
von t von der Gropenordnung |B| bleiben, d.h. fiir die 


(15) [x()| S Const. |B| 
gut fiir alle t; z.B. die eindeutig existierende kleine mit P pertodische Lésung: 
e(é+ P)=2r(8). 


Bleibt noch der Ausnahmefall tibrig, daB zwar (13) mit P periodische Lésungs- 
vektoren 3; (2), ..., 3,() besitzt mit 1<7<n, daB aber fiir alle diese periodischen 


1 Vgl. hierzu die Arbeit [2] des Verfassers. 
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Lésungen 
t+P 


(16) Sie fdt=o0 (eA ff) 


gilt. In diesem Fall lassen sich ohne weitere Voraussetzungen iiber die Funk- 
tionen g;(u, ¢) keine allgemeinen Aussagen machen. 


§ 2. Eine Umformung 


Setzt man 
(17) u(t, A) =ug(t) + Art), 
so ist 
(18) UC, 0) =ato(t) so U(t, Mieed, 24. v(t). 


Wendet man auf die Differenz rechter Hand in der 7-ten Komponentengleichung 
von (9) die Taylor-Entwicklung hinsichtlich 4 an, so kommt 


1 
dg; t d2g; DANE. 
+ gly +E, #) = giltg, f) + Eton) 4 f SHED 4 ayaa 


0 
oder, wenn x,(t) die Komponenten des Vektors x (é) sind, 


“ 0g; ott 
£ phe) 8. (a= 2, Bi (Ug ) oe, 


OUR 


St ei u(¢, A), ¢) al 
+13 et mm (1— A) da. 


(19) 


Wir fiihren nun neben der Matrix at aus (11) noch den ,,Tensor‘ T(u(, A), 2) 


ein, dessen x Komponenten die Matrizen (k, /=1, ..., 2) 
Be /EC2 2 (ued) ae) een 
(20) I; (u(t, a), ¢) ( ate ) jee sid 
sein sollen: 
%, (u(é, a), 2) 
x t t 
(21) T(u(, a), 2) =| *? (u(t, 4), 2) 


%,, (w(t, A), ¢). 
Dann 1a8t sich (9) in der Form 


a2, 2 xr + AiO + f2TOT WGA. 20 1—Aaa 
schreiben. — Hiervon werden bei iter at | kleinem ¢ und f, Lésungen gesucht mit 
(23) |e@|Se bei |B| Sfp. 

Es sei noch angemerkt, da der Tensor 


T (u(é, 0), 4) =T (up), 4), 
die Periode P besitzt: 


(24) T (utp (+ P),t-+ P) =T (up (2), ). 
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a 


Aus dieser Tatsache in Verbindung mit der ersten Beziehung (23) folgt die Exi- 
stenz einer endlichen Konstanten M derart, da fiir das in (22) auftretende 


Integral 
(25) 


gilt. Dazu braucht nur vorausgesetzt zu werden, dal alle rechter Hand in (20) 
auftretenden zweiten Ableitungen in einem Intervall 4; </<¢#,+ P fiir alle Werte 
u(t, A) — 9 (t)| Se des ersten Arguments beschrankt sind. 


nao T (u(é, A), t) x @ (1 — A) ai < M- (max |x (¢)|)? 


§ 3. Der Resonanzfall 


Das adjungierte System (13) besitze einen mit P periodischen Losungsvektor 
3(t), fiir den (14) gilt. Im Gegenteil zur Behauptung des Satzes 1 machen pt 
fiir eine beliebige Lésung r(¢), die der ersten Voraussetzung (23) gentigt, zusatz- 
lich noch die 


(26) Annahme: |r()|<NJ|A| 


fiir alle ¢ mit einer spater zu bestimmenden endlichen Konstanten N. Analog 
zur Formel (33) in [2] erhalt man aus (22) und (13) fiir beliebiges positiv ganzes 
m rit C aus (14) 


3° () [z (t+ m P) —x()|=mC-B+ 


(27) t-+-mP il 
== Ja" (a) fet(e) Tule, A), t) x(t) (1— a) dadt. 


Wegen der Periodizitaét von 3 (¢) gilt fiir alle ¢ eine Beschrankung der Form 
(28) la @| SD. 


Benutzt man (28), (26) und (25), so kann man nach (27) so abschatzen: 


D-2N |B|?+ DM N?|B| m P=m|B||C| 
oder 
(29) 2DN=m|Bp|*- B 
mit 
(30) B=|C| _DMN2P>0. 


Setzt man fiir N eine hinreichend kleine Konstante ein, so ist B positiv. Nun 
kann man zu jedem £ =-0 die Zahl m so gro8 wahlen, daB (29) einen Widerspruch 
enthalt. Damit ist die Annahme (26) in Verbindung mit der ersten Beziehung 
(23) als unméglich erkannt, und es gilt der folgende 


Satz 3. Im Resonanzfall gibt es vier Konstanten 2, M ,D, N, aus denen sich 
vermoge (30) ein positives B bestimmen lapt, derart, da jede Losung x(t) von (9) 


in jedem Intervall der Linge ar 7 P mindestens einmal einen Wert annimmt, so 


dap B| 
(31) |x ()| > Min (e, N V|B]) 
zu Recht besteht. 
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Zusatz. Beschrankt man f durch 


ie eae mit Ninel (vgl. (30)), 


so gilt in jedem Intervall der Lange re P mindestens fiir einen Wert von ¢ 
2 


sogar die Abschatzung | 
(32) JeO|>N VIB. 


§ 4. Der Hauptfall 
Wir bendtigen im folgenden den Hauptsatz iiber implizite Funktionen als 
Hilfssatz 1. Der Vektor 


or Seen eh 
33) Di les 0 (le gecko dag) eit 2 inte 
Up (yy «+5 4y,B) 
besitze in einer gewissen Umgebung der Werte 
WO; dg 0, an, = Op = 0 


stetige erste Ableitungen nach den a, und sei in allen +1 Verdnderlichen stetig; 
auBerdem sel 


(34) (Os cx. 0;.0) =D 
schlieBlich sei die Funktionaldeterminante 


00; (0, . On O) 
Op 


(35) SEO PUA) he ep ree iee 


Dann gibt es zu jedem geniigend kleinem ay) eine Schranke f,, derart, daB zu 
jedem f mit |£|<f, eindeutig ein Loésungsvektor a=a(f) des Gleichungssystems 


(36) v(a, B) =0 


existiert, fiir den |a|<ay ist. Zum Beweis vgl. z.B. [3]. 
Die Anwendung dieses Hilfssatzes geschieht so: Man charakterisiere jeden 
Lésungsvektor x(¢) von (9) durch seine Anfangswerte: 


ay 
a 
G7) (0) estes 
dy 
also 
(38) r=r(a,B;¢). 
Dieser Vektor hat beziiglich ¢ dann und nur dann die Periode P, wenn 
(39) v(a, B) =z (a, 8; P) —x(a,B; 0) =0 
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ist; dies entspricht der Behauptung (36) des Hilfssatzes. Fiir B=0 ist der Null- 
vektor Lésung von (9); also ist, da dann auch a der Nullvektor ist, (34) erfillt. 
Die Ableitungs- und Stetigkeitsvoraussetzungen des Hilfssatzes sind nach be- 
kannten Sdtzen iiber die Abhaingigkeit der Lésungsvektoren von (9) von den 
Anfangswerten a und dem Parameter f gleichfalls gesichert (vgl. z.B. [4]) durch 
folgenden 

Hilfssatz 2. Im Intervall 0<f<P kann man zu vorgegebenem ¢ je ein 
positives ¢, und f, bestimmen derart, dab die Lésungen x (#) von (9) im ganzen 
Intervall der Abschitzung |x (t)|<¢ gentigen, sofern nur |x (0)|S¢, und |£| <A, 
gewahlt wird. 

GemaB (39) muB zu Anwendung des Hilfssatzes 1 noch dafiir gesorgt werden, 


daB sowohl ¢,<a, als auch £,<f, ist. — Differentiation von (9) nach q, liefert, 
da die Vektoren 

a) wyaia 
(40) yt () = He (}= 4,2, .0., 0) 


Lésungen des Gleichungssystems 


(44) Y= Bi) y* 

sind mit der Matrix 

(42) Bp) = (AEM FEOD | GRAM). 
Daher sind die Vektoren 

(43) 9) = ABO OOD bat.) 


Lésungen von (10), da dann (42) in (41) ttbergeht. Die Vektoren (43) sind linear 
unabhangig, weil ihre Determinante fiir ‘=0 die Einheitsdeterminante ist (vgl. 
hierzu [2], Satz 4). Da wir hier im Hauptfall die Voraussetzung machen, daB 
(13) und damit auch (10) (vgl. [2], Satz 8) keine mit P periodischen Lésungs- 
vektoren besitzt, so ist (vgl. [2], Satz 10) 


(44) Det] MIE” | 40, 


Fin Blick auf (39) lehrt, daB dies gerade die Bedingung (35) des Hilfssatzes ist. 
Da alle Bedingungen des Hilfssatzes als giiltig nachgewiesen sind, ist seine Be- 
hauptung, bei uns also (39), bewiesen. 


Es gilt daher 


Satz 4. Hat (13) keinen mit P periodischen Loésungsvektor 3(t), so gibt es zu 
vorgegebenem & Schranken e, und B, derart, dab zu jedem B mit |B| SB, ein Anfangs- 
vektor x(0) mit |x(0)|<e, existiert, so daB die Lésung x(t) mit diesem Anfangs- 
wert x(0) periodisch ist und der gleichméfigen Abschatzung |x (t)| Se geniigt. 
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Dies ist ein Teil des Satzes 2. Um Satz 2 vollstandig zu beweisen, bendtigen 
wir noch den 


Satz 5. Es existiert eine endliche Konstante E derart, dap bei geniigend kleinem 
|| fiir diese mut P periodische Lésung x(t) von (9) gilt 


(45) : |S E- |B]. 


Beweis. Sind 3, (¢), ..., 3, (¢) linear unabhangige Lésungen von (13), so erhalt 
man unter Ausnutzung der Periodizitat von y(t) analog der Gleichung [2], (38) 
Uap —l te 10 

Zt 


t+ P 


Feo} £* (x) T(u(z,A), 2) £(t) (1— A) dade 


(46) 


(vgl. dazu auch Gl. (27)). Die Determinante der Koeffizientenmatrix linker 
Hand ist nach der Voraussetzung des Hauptfalles von Null verschieden (vgl. (44)). 
Daher kann man das lineare Gleichungssystem (46) fiir die linker Hand stehenden 
Komponenten %,(é),..., %,(¢) von x(¢) nach der Cramerschen Auflésungsformel 
auflodsen. Bezeichnet bei festgehaltenem ¢ 


7 == Mak se ersiany FE) ig 
(47) x = Max (|%,(t T)| (v ' n) 


so findet man auf diese Weise eine Abschatzung der Form 


(48) 


- (max |x (¢)|)? 


mit zwei endlichen Konstanten A und L, wobei (25) beachtet wurde. 
Da 


ist, folgt aus (48) erst recht 


max |x (¢)| S/n K 


\nL- (max |x(6)|)?, 
also a 
max |x (¢)| [1 — mL - (max |x ()|)] S/n K| BI. 


Auf Grund von ane 4 diirfen wir OL aay daB die eckige Klammer positiv 


ist. Fires Ae z.B. wird sie > . Daraus ergibt sich sofort die Abschatzung 


(45) mit E=2 nk. 

Statt die Cramersche Auflésungsformel anzuwenden, kann man natiirlich auch 
(46) durch linksseitige Multiplikation mit der zu [an (¢+ P) —3, (¢)] inversen 
Matrix nach r(¢) auflésen und dann mit Hilfe der Schwarzschen Ungleichung 
|x (#)| abschatzen. Auch so gelangt man zu (45). 


II. Teil einer bei der Naturwissenschaftlich-Philosophischen Fakultat der Tech- 
nischen Hochschule zu Braunschweig eingereichten Dissertation. Referenten: Prof. 
Dr. R. Ieriscu und Prof. Dr. H.-J. Kanotp. Auch Herrn Dr. H. ELTERMANN bin ich 
fiir wertvolle Ratschlage zu groBem Dank verpflichtet. 
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Karim: Resonanz bei nichtlinearen Differentialsystemen 
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§ 1. Introduction 


The subject of this paper! is the development of a numerical procedure for 
finding an approximate solution of the Goursat Problem: 


(1.1) Uyy =] (%, Y, U, Uz, Uy) , 
(1.2) u(x,0)=ao(x), “(0,y)=t(y), o(0)=7(0), 
(1.3) OEao r= ds OS Ws On. 


We shall denote this problem by 2-HP. 
The usual Runge-Kutta procedure, which we denote by R-K 1, applies to 
the problem of finding a function y(x) satisfying 


dy 


dx =f (%, 9), WO = Vee Ore LS a, 


(1.4) 


The procedure developed here, denoted by R-K 2, is analogous to R-K 1 
both in mode of application and order of accuracy with respect to Taylor series 
expansions. Just as the usual Runge-Kutta procedure is a refinement of the 


1 The research presented here was supported by the Office of Scientific Research, 
Air Research and Development Command, under contract No. 49(638)-192 and by 
the United States Naval Ordnance Laboratory. 
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Euler-Cauchy polygon method for (1.4), so also R-K 2 is a refinement of Draz’s [1)} 
analogue for 2-HP of the Euler-Cauchy polygon method. Again, just as the 
validity of R-K 1, in the sense of uniform convergence to a true solution, depends 
on the Peano existence theorem for (1.4), so also the validity of R-K 2 depends 
on an analogue for 2-HP of PEANo’s theorem, given in [3]. 

Concerning regularity, it will be assumed that f is three times continuously 
differentiable with respect to any combination of its arguments, that / and its 
first derivatives are bounded, and that o and 1 are five times continuously 
differentiable throughout any domains in which they are used. Naturally, if 
a lower order of accuracy than that given here is acceptable, then these regularity 
assumptions may be relaxed. While these assumptions on o and t are necessary 
in order to assure the stated accuracy, only tabulated values of o, 0’, t, and t’ 
at suitable points are needed in order to carry out the actual numerical com- 
putations. 

It is suggestive to consider the system of integral equations equivalent to 
Ci) and(422): 


u(x, y) =a(x) + 1(y) — (0) +S JHE"), En), 9m) dé dn, 


es )) 0.2) a: f f(x, u(x, 1), B(% 0), 9(% 1) an, 


g(x,y) ='(y) + ie y, (Ey) bE. 9), g(E, 9) dE. 


The procedure to be described can be viewed as a technique of stepwise 
approximate integration of these equations. 

Thus, the overall numerical solution of the problem in a rectangle R: OS xSa, 
O<yS4, is to be carried out in a stepwise manner over a rectangular mesh on R. 
The object of R-K 2 is to provide a method for obtaining approximate values 
of u, wu, and uw, at (x»+h, Yot+qh) given the values at (%9, Vo), (% +, Yo) and 
(%9, ¥o + gh) where these points are the corners of a rectangle of the mesh. 
In terms of the Taylor expansions of u, =u, and g=u, about (xp, Vo), these 
values are to be accurate through terms of order fh‘, h® and h’, respectively. 
This implies tacitly that w is bounded. 

The solution of the differential equation (1.1) on the mesh subrectangle 
determined by (%, yg) and (%)+h, yo+@h) is a “miniature” problem of type 
2-HP. For the derivation of R-K 2, one assumes that the values of w, u, and u, 
at (%, Vo), (X» +A, Yo) and (%, Yp+h), as given from computations on the 
preceding subrectangles, are exactly the true values for the solution of 2-HP 
on R. The functions s(«) and ¢(x) defined by 


(1.6) S(x) = u(x, Yo), — t(y) = U(X, 9) 


serve as the initial conditions for the miniature problem; they have the same 
differentiability properties as the true solution uv, so u,=p=s’ and u = 7 
The original initial conditions given by o and 1 are involved only inaplicitly in 
the formulas. However, for subrectangles along the initial lines (those adjacent 


tOp%=—0) and y=0), explicit use will have to be made of o and t (cf. eqns. (3.8) 
and particularly §11). 


| 


| 
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The basic idea of the derivation is then the same as that for R-K 1: An 
explicit computational procedure is outlined in which there are several free 
parameters; here these are the elements of certain diagonal matrices (cf. equa- 


tions (3.13)) as well as quantities permitting consideration of only part of the 


subrectangle with sides 4 and ph (cf. §3). The values of wu, p and g at (x%)+h, 


-Vo-+qh) are expressed as power series in # in terms of these parameters. Then, 
_ by matching the coefficients of the lower order powers of / in this expression 


and in the Taylor expansions of u, # and g about (x9, Vo), a system of equations 
is derived with the parameters as unknowns (cf. equations (R1)—(R 11) in 
§§ 6—8). A solution of the system, used in the computational procedure, will 
then yield values for u, p and gat (%)+A/, vp + wh) of a prescribed order of accuracy. 


§ 2. Taylor Expansions 


Under the assumed regularity of /, o and t there is a unique solution w(x, y) 
of 2-HP. With this solution, let 


F(x, y) = (4%, 4(%, ), Ue (%, Vy), My (%, ¥))- 


Then F,, F,, F,,, etc. denote the partial derivatives of F. Note that, from the 
differential equation, f,=wu,,,, for example. 
If one writes the Taylor expansion about (x», Vo) for w(%+h, Voth) and 


makes use of the last remark and the definitions of s, ¢ and /, then one finds 


Au =Uu(% +h, Vo + ph) — U(X, Vo) 


bi 3 Lad h wt h* 
a=ihis guile pie eS sip ee st 
2 3! 4! 


, gy? h ’” gp? h> yr gh 
(2.1) (Pa ee ep StL alee OO 
1 1 ~ 
gh? f+ 9 h|—E4 5 9B) + 
+ 9M Et 4 pF, + 5 PB] + O08). 


In this, the various coefficients of the powers of h are to be evaluated at (%, Yo); 
in particular, the / in ph? means F(x, Yo). 
Similarly, for the function p=w,, since p, (x, vy) =F (x, y) and p(x, Vo) =s'(x), 


Ap=f(% +h, vo + ph) — P(%o, Vo) 


(2.2) =hs +43!" | “ s+ ohft@gh? 


E+ +95 ai 


tpl Et Sp Ey t 2 oF] + OM). 


In the same manner, one finds 


3 
Ag=phi'+ 28" yn 4 OM 04 nj +n E+ 95] 4 


(2.3) 1 4 1 
EY W| 2 Et Ens Fy+ 2 ye E,| + O(h'). 
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§ 3. Definitions of the Matrix Quantities to be Used 
As in the one dimensional case, it is not very practical for numerical work 
to use these Taylor expansions directly to determine the increments Au, Ap 
and Ag. Rather, the various terms will be matched by the procedure below 
to within the desired order of accuracy. For this purpose, a matrix notation 
will be employed. 


Let 

OY a 

Z 1 Ox 
, ‘ 0 se 10 of 0 

oe 

Gl) U=("], SU)=| 0) |, F=| 0 OF, 0 0 
p 7(U) Bie tang ler 
q {(U) aS 9 
0 0 0 0 = 
x 


in which /(U) has the value /(x, y, vu, p, g). Then the differential equation (1.1) 
can be written 


(3.2) QU=(V). 


It is to be noted that in general the first two components of vectors and 
diagonal matrices are thus associated with the independent variables x and y, 
while the latter three components are associated with the dependent variables 
u, p and q, respectively. 


The initial values at (%9, Vo) are expressed by the “‘initial vector’’ 


The initial conditions along the boundaries x= x» and y= yy of the subrectangle 
under consideration do not need to be given in continuous form; indeed, because 
the stepwise numerical process gives only values at the mesh nodes, only finite 
differences are known. The vector expressing the conditions along *=%,) and 
y=, is thus taken to be 


0 
@) 
(3.4) B= Ast 
As’ 
At 
where 
As=s(%» +h) — s(x»), At=t(yy+ ph) —t(y), 


As’ =s'(% +h) — 8'(%), At =t(yyo + ph) —#'(y). 
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There are three things which should be noted here which point out the modi- 
fications of R-K 14 needed for 2-HP, and in particular make plausible the choice 
of the increment vectors H? below. First, as suggested by the fact that there 
are three unknown functions, viz u, p and q, three cases will be distinguished, 
one for each of u, p and q. That is, it will in general be necessary to match 
separately the Taylor series expansions of Au, Ap and Ag. (However, see the 
comment at the end of §9 concerning the possibility of other solutions (not 
now known) for parameters satisfying the requirements (R 1)—(R 11) of 
§§ 6, 7, 8.) Matrices and their elements used for a particular case will be dis- 
tinguished by a superscript u, # or q; the superscript w will stand for any one 
of wu, and g. For notational convenience, the superscript will often be suppressed 
if no confusion can result. In particular, if a matrix has a superscript w, then 
each of its elements is assumed also to be so indexed, even though the w is 
suppressed on the elements. 

Second, because in the Taylor expansions for Au, 4p and Ag the terms 
involving second derivatives of F are three in number, enough freedom is needed 
to match all three terms. (In the Taylor expansion of y (x) used in the derivation 
of R-K 1 for y’=f(x, y(x)), there is only one term involving the second derivative 
of F=f(x, y(x)).) There will be three equations for matching these coefficients, 
viz (R4), (R5@), (R6q@) of §§6—8. For this reason, entities will be intro- 
duced with a subscript A where A=1,2,3. The greater degree of freedom 
thus afforded will permit satisfying the mentioned requirements. 

Third, if R-K 1 is regarded as a procedure for approximate integration of 
the integral equation equivalent to the differential equation (1.4), then the 
increment coefficient / used in R-K 1 is seen to be a measure of one subinterval 
of the domain of integration. From the integral equations (1.5) it is clear then 
that the corresponding increment coefficients to be used for u, # and g in R-K 2 
are ph, wh and h, respectively. These appear as the last three components of 
H? below. 


Thus increment matrices H? are defined by first setting 
(3:5) 2 =P and ho = gih (== 16, PD, Os A == 3) 
where ef and g? are parameters to be determined. Then one defines 
hy 
Pail, 
(3.6) H? = diag| 9,5 
Pal, 
h, 
(Note that the w has been suppressed.) 


This definition of the H? then implies integration over only a portion of 
the subrectangle S determined by (%p, Vo) and (%)+%, Voth); that is, sub- 


rectangles S? determined by (x9, Yo) and (% +7, Yot+ p7h%) are also considered 
(cf. §9). By requiring 


(3.7) Cray it andi 4 c7e7, Sits, 
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one is assured that each S® is contained in S (and is nondegenerate). For these 
smaller subrectangles S%, further initial vectors corresponding to B will be 
needed. They are, to begin with, considered as a simple “cutting down” of B 
to S® by means of linear interpolation. Thus s(%»+/7) —s(%»), which might 
be thought of as A?s, is replaced by g?4s. By the mean value theorem of 
differential calculus, this can then be replaced by g?-/5’, that is, by h?s’, where 
3’ is a suitable ‘‘intermediate value” of s’(x). It turns out, in the derivation, 
to be necessary to define also certain other intermediate values §’, 5”, ¥ prs 
(See the comments following (5.18) and at the end of § 6.) The definitions for 
all of these are given by: 


hs’ =s(% +h) — s(%p), 2h 

hs" =s'(%) +h) — s'(%), 25" = s'(%) + Sh) — s’(%), 
ght =tyy+—h)—t(y), Bpht =t(yo+Fph) —t(y), 
phi’ =t(yyo+ ph) — tly), Bpht”=l(yo+ $ph) — ty). 


(3.8) 


The numerical determination of s’, 5”, /’ and ¢” to within the required order 
of accuracy so as to be able to prove the convergence of the approximate solutions 
is one of the more delicate parts of R-K 2. It will be discussed in detail in § 11. 


Using the above, one defines the B? as follows: 


6) 
6) 
Bh =| WS + that’ |, 
has 
(3.9) mi 
O O 
O O 
Be=| Ast giiht |, Bl=| ns’ + eqef At 
gh As’ ngs” 
gh hit!" ele At’ 


These are the best forms of the definition for numerical purposes, but for the 
derivation to follow it is more convenient to replace 4s, At, As’ and A?’ via 
(3.8) so as to write these simultaneously : 
O 
O 
(3.10) BS | Wes nay 
pee” 
99 het” 
in which ~ means — or ~ as given by (3.9) according to which case, u, p or g 
is under consideration. 
2 erences such as s’— s’(x») will be important in the sequel. For brevity, 
ets,s ,?,¢° stand for s'(x), s’"(%9)) t'(yo), Vo) respectively. Then, by two 
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term Taylor expansions, it follows from (3.8) that 


S’— s' =% hs" + OW), g’— s' =8,2hs" + O22), 
gh Sei St O(h?), S72 si= 02h gee ae O(h?), 
f—t =O, pht"+ OR), t—f=8,2pht"’+ On), 


~ 


fh =O, pher+ O(h?), Mie t/’= O,2ph oe O(h?), 


(3.11) 


where 0<%,, O;< 1. This last simple statement about #; and 8, is not specific 
enough for later purposes. It will, however, be enough to use the known fact 
(see, e.g. Prasap [4]) 

(3.12) 0,= 34+ Oh) and 0,=¢+ OA). 


a) 


Besides the parameters ef and g? so far introduced, further free parameters 
are introduced by the matrices (defined for w=u, p,q; A=1, 2, 3): 
VP = diag (Vi1,Vi.2, Vis» Via Vis), 
(3.13) Wy = diag (Wasa, Wise, Wiss, Miia, Mais) (t= 1,2), 
Aj, = diag (Aj;1,Agi2,Ari3,Aria,4as5) (¢=1,2,3). 


Here, 7;; means the 7" component of the vector or diagonal matrix T,. 


§ 4. The Basic Equations for the Runge-Kutta Procedure for 2-HP 

As mentioned, three cases are to be distinguished, one for each of u, p and gq, 
in the following procedure. In each case three sets of computations are to be 
carried out, one for each of A=1, 2 and 3 and used in equation (4.3). (However, 
in the actual numerical work, the particular system of parameters found in 
§ 10 permits one set of computations for each of the cases for # and gq to be 
eliminated. In the following derivation this will be ignored.) 

The equations below setting forth the computational scheme are very similar 
to those for R-K 1; indeed they are almost identical if the equations for R-K 1 
are written in the analogous matrix notation. Besides the indices w and 4, 
the principal alteration is the introduction of the B and B? matrices which 
are not present for R-K 1. The equations are 


K?, =H? ®(U,), UA =U) + (Ki + Bi), 
(4.1) Kj, = HP @(Un), Uji2z= Up +E ai (K7; + Ba), 
Kj, = H? ®(U72). 
To define the weighted average corresponding to that in R-K 1, first let 


h 


ph 

(4.2) KG,=|"Kii3 |) 
KGa 
K4i5 
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and then define 


(4.3) Ae | eA Ks) +B. 


The approximate increments in u, p and g are then taken to be the last 
three components of A, respectively: 


U(%q +h, Vo + ph) = (xX, Vo) +A, 
(4.4) b(%p +h, Vo + h) = (%>Vo) + 4a, 
q(% +h, ot Gh) = q (Xo, Vo) +A. 


§ 5. Computations Preparatory to Matching the Taylor Expansions 


It is the purpose of this section to derive detailed expressions for the elements 
in each K,, and, further, to examine the differential operators which are employed. 


The superscript w will be suppressed throughout the section. 


To begin with, recall the general Taylor expansion for a function of five 
independent variables: If the differential operator D is defined by 


then2 oes 
co 


(5.1) Fat 1s %2+ 72, Wg + 1%, %q + 74, %s +75) = ee ace Oo, Hy, %4, %5) « 


n=0 
It should be noted that the values of /, F, s,¢ and their derivatives in the 
following equations are those at (%», ¥p) unless otherwise indicated. 
From the definitions in §§ 3 and 4, one has 


h, 
Pal 
(5.2) K,, = 1, ®(U,) = Pp hi, f 
Pil, t 


h,f 
and 


Xo t Vash, 
Yo t+ Va Gah, 
653) U,=U+V(Ky, + B,) =| u%+ Visit ths’ + YP h,t’) 
Pot Va a(Galal + 8") 
Qo t+Vi shat + o,h,t") 


Let 
(5.4) Vg= 0, = Ur os Aes 
and let 
(5.5) Dj, =v, DF 


27 j 
f necessary the series can be replaced by a truncated series with remainder 


term. Because of the regularity assumptions on f, e 
i , enough 
the analysis below is not altered, ‘ : Cee jme cre het Ee awe 
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where a a hte 
Diy ——+ (wmyf+s'4 gail’) s 
(5.6) “4 Sf M3 
1p (a) A ts) 
(pa f S ) Oc, (f Q,t oh 


Then in the Taylor expansion of f(U,,), 4, D,, plays the part of D in (5.1). Thus 


TU ee Se D2 7, 


From this n=0 
h, 


Pal, 
- De hi 
pitt hDaft— Dutt Duit--| 


(5.7) KK,» ae HH, ®(U,,;) Ty 


’ 


pee 
oe oe 


h h® 
hlf+hDat 4 Diet : (Di jit | 
so that 


2 
U,..= + 2, Wis (Has + B,) 


xq + (Wars + Wyo) Ia 
Yo t+ (Wire + Wioe) Pala 7 
(5.8) =] % + Wars Kars + Wios Kias + (Wars + Mras) a(S’ + Gal’) 
Po + Wira Kara + Waa Kana t+ Mara + Mia) ns" 
Go + Wars Kars + Mies Kaos + (Wars + Mies) Galat” 


Let 

(5.9) Was tWiay=es 4,555 2=1,2,3), 
and let 

(5.10) Daw, DF: 


Then, with (5.8) in mind, one can write, using (5.2) and (5.7), 

yh, + 10, Pa hy sa + Mars Kars + Wars Kins + wiles ‘+ yt) 8 = 
+ [Wyia Kara + Waos Kao, + Oy 8"] ae a 

+ [Wars Kars t+ Mies Kies tm gily ca 52 

= hy Din + Wins (Kina — galtsf) G+ 


0 5 
+ Wryoa(Kae4 — Gila) ae + Wyos (Kis — Mf) 


Oa; 
h he a 
=I, Dj + Wars al Daf + Dist x Dif Aes + 
2 hy n9 hi 3 (4) 
+ Whoa Pall Dia f+ if 3 | Dif oe 


h hi 
+ Wins li| Dar f + Di if- a Din | f--| 5 
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It is this entire operator which plays the part of D in the general Taylor expansion 
(5.1) of f(U,.). Thus 


(U2) =f+hDiof+ 


ae 14 Dat + Wena Dar f- fa + Wins Daal fel + 00), 
SO a 
DP, hy, 
(5.12) Kj aU ee 
Pp h, f (U;») 
hy, f (Uy) 


It will be seen in equations (6.1) and (7.1) that the operators D,, and De 
are to be compared with derivatives of F. In order to carry out the comparison, 


let {; denote Sc Ha Gy, Wig On) (== 4 palees pcm aee 


ou op og 
i is ae ee Wilders 
The values of interest for the functions in the foregoing computations are those 
at (%), Vo). At the point (xp, Yo), because of (1.6) and (1.1), one may write 


(5.13) Fe=At+s'fpts" fat fts- 
Similarly, at (9, Vo); 

een BR=ftl htt ts 
and 


LS face S Asses Baa tise | 
+ 2(8' fis +s" tat fhe $s" Isa 8 Fiestas thas) + 
Sifa eS laches 
Fey =hets'U fgg +8" faa tt" fisst 
Uhsts'fastthatsfeatt fs tiifest 
F (Sf ES") fgg + (80 +f) fg + (87 2" +P) fas + 
fis t+Litat Fyfs, 
Fy = hee t b' fs5 4 Pigs a 
+20 festtheatt’ fos tl tiaa ttt” fos +t" thas) + 
Pesce hy gar’ fee 


(Note particularly the terms in fj, f, and f;; they are the root cause of 


the definitions of the s’.¢s".#/~ Their presence here means that, for later 


purposes, the §’,#’,§”, ¢”” in (5.15) below may not be taken as s’, #’, s’’, t” 
respectively; cf. equations (6.2), (6.3) and (6.5).) 

From the form of (5.13) and (5.14) and the definitions of D,, and D,, it is 
seen that f+, 4, and D# f should be compared. Indeed 
E+ pak, — Df 


=[—aht— 6 —8)— 9, =), = © — 39 f= SP) oa. 


3 


(5.45) 
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Distinguishing now the meaning of ~ for the cases for w, p and g according to 
(3.9), one has from (3.11) and (3.12) (still suppressing the superscript) in the 
case for u: 


(5.16) Be Git, — Dz ft 

3 —[(piinf thks" tigi pht’) fytths' f,+4qiphe” fe] +O); 
in the case for p: 
(5.17) 1S te Dee 

Bt We 

=—[(giiitshs’+ sei Pht’) fet ths’ fytig mht” fs] + OR); 

and in the case for g: 
(5.18) LE a Ap Rey Df 

=—[(@mttghs +e vht) fpt+shs’” fatto, phi” fs] + O(n). 


To have the coefficients $} and } appear as shown was the object of introducing 


the quantities $’, s”’ , t’ and ¢” into the B’. The reason for having these particular 
coefficients is nord at the end of § 6. 


The above equations (5.16), (5.17) and (5.18) each represent three equations, 
one for each of A=1, 2,3. Later, we shall wish to multiply the equations (5.16) 
each by a certain «, and add, and to do the same for (5.17) and (5.18). 


Note further that 
(5.19) ote 


because the coefficients of the d/¢«; in D¥ are constants. Thus, for the second 
order operators, the comparison is between D¥” and F,,+29,F,,+ 95 Pe leh 
equation (6.1)). One finds that 


DF? f= hit 2Qrhet Vifeet+ 

+ (Pilaf +S + git) fas + Pat +8") feat f+ Gal)? hos + 

+2[ (putts + Pat ieee ies Nae (f+o.0") fis + 

+ P(Galat+ §' Pat ) fas + pulpal +8”) faut Gall + Gat") fas + 

+ (pilaf +8" + gat ) (Pat +8") faa + (Qilgf+s’ + qt ) (f+ gat’ ’) tes + 
+ (gaf +8") f+ gat’ PNA 


For the comparison just mentioned it will be sufficient in the sequel for terms 
merely to agree within order 4. Thus, recalling that s’=s’+ O(h), etc., one has, 
for example, for the coefficients of fg, that 


(s’ + git’)? = (palgt+s'+ gat y!)2 O(A), 


and hence 
vet 2Q,F, yt GiFy — DFT 


(5.20) = (s+ 2y,f+ 950") f+ (s+ 2g, F, + 07%) fa + 
+ (F,+2m,F,+ git’) fs+ O(n). 
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As with (5.16), this equation stands for three equations: A=1, 2, 3. Ot the 
coefficients of f;, f, and f; here, those parts involving s’, 7’, s’”, t” are similar 
to those obtained in (5.16—18). It is this fact which, through the definitions 
of the B®, allows the complete matching of terms in the next sections, as in the 


computation of (6.5). 


§ 6. Requirements on the Parameters in the Case for u 

In this and the succeeding two sections, requirements on the parameters 
analogous to those in R-K 1 will be derived. Since in this section only one case 
(viz u) is considered, the identifying superscript will usually be suppressed. 

From the equations (4.1), (4.2) and (4.3) defining the computational procedure, 
one has 
A= Ds Ajis Kiss) + Bs. 
mer 
In A, the quantities KZ; 3 and Bs; are to be replaced by their equivalents as 
given by (5.2), (5.7) and (5.12) (and (5.11)) and the definition (3.4). Here, the 
quantities A s and A¢ involved in B.; should be replaced by their Taylor expansions. 
The terms of 4, are to be grouped according to powers of / and compared (by 
subtraction) with the corresponding terms in the Taylor series expansion (2.1) 
of du. The terms of (2.1) involving s and ¢ are exactly those terms in the Taylor 
expansions of 4s and At which appear in 4, because of B,. Hence, all the 
terms in s and ¢ cancel when A, is subtracted from Au. The remaining terms 
of 4, involve the formulas for the K,; 3. In these it is to be recalled that 


i= 4, hy = gyh, 


so that m and h can often be factored out of a given term. Thus, letting ¢, 
denote the absolute difference between Au and 45, one finds 


é3=|du—A,| 
=19 h*7 1D Arsen a 


r 3 
(6.1) + Phila, + PE, — 2 (Arese8iDah + Aj 33 e188 Draf)| it 


A 3 
ph TOL. LR 21 (Anas 8%: Dis f+ Arsseagh BDia fh) — 


3 
ee e4 Bi (Wie PiPart- fea t+Mios Dart: fs)|| + O(n’) . 


The object is now to find requirements on the parameters such that, if the 
requirements are satisfied, then ¢, is of order h5. This will be done by demand- 
ing that the coefficients of the lower powers of i in (6.1) be zero. The coefficients 
will be treated in turn. First, to eliminate the h2-term requires 


(R 1 u) 2 Anse ot 
3 
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For the coefficient of h®, recall that —Dia= ’ Dis Di Thus tD7\ isea 
A cy 


factor of each term of the summation in the coefficient of 43. The coefficient 
of Dj is taken as the definition of an “intermediate” parameter O43: 


(A;) Hag = & 3 (Ajo3%at+ Ags W). 
Now require 
3 
(R2u) Diss; 
ail 
3 
(R 3 u) Dos =F. 


Note that the coefficient of F, involves » in the 3-term of (6.1), while in (5.16) 
it involves y,. This is the reason for the e, which appear in (R3u). Multiply- 
ing the A equation (5.16) by «3, adding the resulting three equations for 
A=1, 2,3, and making use of (R 2u) and (R 3 u), one thus sees that the coefficient 
of mh? in (6.1) is given by 


3 
BE +495 —YiajsDt 


(6.2) =—|(Des pital beaters b Spr ht’) fy + 
+ EES Ath Egthe” fs] + OU). 


The right member of (6.2) has / as a factor; since the member is already part 
of the /?-term, the coefficient of # thus found here can be moved into the coef- 
ficient of h* in (6.1). The O(h?) term in (6.2), when multipled by 43, can be 
incorporated into the fifth-order term of € 5. 


Hence, the coefficient of wh* in (6.1), including the part from the /*-term, is 


eee oe Jp asa) ~ 23 ea 84 (Aj2g%%4+ Aas 0%) D3? fl — 
(6.3) =a 2 Aj33 2 81(Wios Pa Dart fat Wios Darl - fs) — 


_|Sasnaltastier|atentaer' i), 


Let 

(B,) Bas = 2 84(Ajo3 04+ Arga W%) 
and require 

(R 4 u) 2B; 2 - 

(R 5 u) 2D Bist = 

(R6u) DPrst=F. 


Note that the coefficients of &, and F,, in (6.3) involve y, while in (5.20) they 
involve g,. This is the reason for the e, and e; which appear in (R 5u) and 
(R 6u). Multiplying the A equation (5.20) by £,.,, adding the three equations 
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for A=1, 2,3 and making use of (R 4u), (R5u) and (R 6u), one has for the 
part of (6.3) in curly braces 


Eat 3 bry4 oP Ey, — LBs Dit 


af 

6 

(6.4) =(is’+loftepiyht es’ +rekh ter iat 
1 (LR +19 Fk +3970") f+ Of). 


KX 


Since the @(h) term in (6.4) can be incorporated into the fifth-order term of 
€3, (6.3) becomes 


(Rei Uesmsllh+teb+3P BH) + GE + EPH) fs — 


(6.5) 
ai pa Ajs3 €,85(Wioa pPDirt fa t+ Wios Dail -fs)- 
2 


(Note the elimination of the terms ins jt 7s' andy.) 
Consider the coefficients of /, in (6.5). Let 


(C,) AES ae e; £4 Ajss Whoa, 
and require 

(R 7 u) 2 Yas 

(R 8 u) Las a se 


Then from these requirements and from equations (5.16) it follows (in the same 
manner as for equation (6.2)) that 


(EE +S0 8 — Sadi i f= Oh). 


Thus, this difference can be incorporated into the fifth-order term of «3. 
Consider the coefficients of f; in (6.5). Let 


(D,) Ors = &% Si Ariss Was 
and require 
(R9Ou) LOs=s 
a 
(R 10 u) Dd 013% =¢4- 
v7) 


Then, again, by adding equations (5.16) each multiplied by the appropriate 6, 5, 
one finds 


(bE +19 F — DdisDt ff, = O(h) 
which can be incorporated into the fifth-order term of ¢,. 
Finally, so that the coefficient of f, in (6.5) becomes zero, require 
(R11 u) Dias Sa=F- 
A 


If all the above requirements (R 1u)—(R 11u) are satisfied, then it is seen 
that ¢,= O(h5), as desired. 
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It is worth noting that the particular coefficients 4 appearing in the right 


member of (6.2) were originally set up by the use of §’, §”’, ¢’ and t’’.in the defini- 
tion of BY. This was, of course, the reason for so ve BY. In the cases for 
fp and gq still to be considered, the same reason will appear for the definitions 
of BG and B34. 


§ 7. Requirements on the Parameters in the Case for p 


The derivation of the requirements in the case for # is, with two exceptions, 
exactly the same as in the case for u. The first exception is that the power of 
h is lower by one in the various terms to be matched in A, and the Taylor 
expansion of 4. The second exception is that equations (5.17) are used instead 
of (5.16). Because the cases are so much alike, most of the explanatory comments 
will be omitted here. 


From (4.1), (4.2) and (4.3), 


Substituting into this the formulas ion § 5 for K,;, and the definition of B,, 
one finds 


és=—|Ap—A,| 
=|pri| — DAvia eos ae 
(7A) +MlF +405 Se (Area ta8h Dart + Ansa %s 8% Dzaf)| + 
| 
I 


pW lh EA ZO Le EG Lyy 2 (Aiea ei Dart t Aisa %s 8a Dial) — 


— Y Anas 84 Wizs Pi Dart fa + Wias Dart -fs)|| + O04). 
a 
To make the coefficient of g 4? vanish, require 

(R 1p) Dd Araahi=1- 

A=1, 2,3 

i=1,2,3 
Next, set 
(A,) O44 = €,82(Agoa Va + Ansa), 
and require 
(R 2p) ama =1, 

A 

(R 3p) >) 4 =F. 


From equations (5.17), the coefficient of mh? in (7.1) becomes then 


F,+4¢9F,—d%4D7 f 
2 


(7.2) == |(Sonegahif+ehs” +$-k9tht") ft 


+ hh +d bpthe'” fo] + (0). 


4* 
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Note that / is a factor of the right-hand side. Thus the coefficient of ph in (7.1), 
with the additional terms from (7.2) included, is now 


{Ped eee, ao" Ey — 3 bey gh(Araa 0h + Ansa’) DEH — 


(7.3) — » Aja €,23(Wro4 Pj Diait-fs+ WMios Dirt -fs) — 
Z 
—|(Deagisif ths s’ +6 PU \is+% eo os +4 g*t”' Fe]. 
Let 
(B,) Bra =F 8i(Aroa Vat Anza wi) , 
and set 
(R 4p) Dia 
(R 5 p) 2 Bsa =o 
(R 6p) 2D Bisi=s- 


From (5.20), the part of (7.3) in curly braces becomes 
yee I Oya EP ale DBra Diet 
(7.4) =(s"+ 329i + eee it Gs" +39 Fk, +eeh) tat 
t BE +39 B+ eet’) fs + OA). 


Then, from this, (7.3) becomes (since the @(h) term can be incorporated into 
the fourth-order term of és) 


(Pei — Lea gisllt GPR + teh) + bE +495) f— 
— 2 Arse 4,85 (Wasa Pi Diart-hh+Mies Dyyf + fs) - 


(7.5) 


Consider the coefficients of f,. Let 


(C,) M4a= GA, aa WVio4 04 
and require 

(R 7p) LV4= 2 

(RK 8p) eZ ea 


From these conditions and from equations (5.17) it follows that 


o(tF 4 cn oes Lia D¥ f\ r= O(h). 


Thus, this difference can be incorporated into the fourth-order term of oye 


Consider the coefficients of f; in (7.5). Let 


(D,) Oy a= 184 Ags Wie %; 
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and require 


(R 9p) 2042: 


A 


(R 10 p) 2 O14 =F: 
2 


From these conditions and from equations (5.17) it follows that 


(2E +29 5 — Dibra DEA) f= Ol). 


This difference can be incorporated into the fourth-order term of e4. 
Finally, so that the coefficient of f, in (7.5) becomes zero, require 


(R11 p) Dia & 8a = 2- 


53 


If the parameters satisfy the above requirements, then ¢,= @(h*) as desired. 


§ 8. Requirements on the Parameters in the Case for q 


The procedure for this case is the -same as before. For this reason only the 


resulting requirements are given here: 


DieAge e415 
3 


=P, 
i=1,2,3 
2 
Oy 6 = 81(Agos Ya t+ Aras Ws)» 
9 
> 01.5 — 2> 
A 
Dea Cis 


2 
Bis =281(Aro5%at Aras wi) , 
2 Bas 
2 Bis A=, 
2 Ba.s ej 


|e 


, 


I 
bole 


3 

Vas — ©a8a AjssWieata, 
D> Vas 
A 


wag ite 
a 
a 


ab 
Q> 


| 


3 
64.5 = 81 4Aiss Wes va, 
» Ons 
2 
D>; O5sa= 
2 


<1 
D5 aka s- 
a 


1 
6? 


| 


dolH 


, 


(Zao +ph) 
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§ 9. Further Requirements 
Besides the conditions on the parameters imposed in the preceding sections 
there are also the natural requirements (see (3.7)) 


O< et, Osc ee ey ate 


It is also reasonable to impose simple requirements on the e, and g, which 
yield some kind of symmetry among the subrectangles with dimensions h, by 
,h,. Each of these subrectangles has one corner at (%», Yo) and the diagonally 
opposite corner at (%)+/,, Yot@al,). The symmetry requirements which seem 
most natural are that for each case, u, p, and gq, (1) the first rectangle (A=1) 
have the corner (%)+/,, Vo +G/4) on the diagonal of the basic subrectangle 
(see Figure 1), and, (2) the 
corners (%9+ Ag, Vo + Po Mo) 
and (%+/5, Vot@3M3) of the 
second and third subrectangles 
be symmetrically placed (as 
dare Wr ee es ah : in the Figure) with respect 
to this diagonal. This symme- 
try is expressed by the pro- 
portions 


Pyhy 


Pp hg 


aR rh 2st’ P2 Nz Pahs 


wo h, hs h, 


From these it follows that 
€y=1, €s=1/e, and gg=—e, Bo, 
with ¢,, g, and gz still free. 


(Zo, Yo) hy hg hy (Zy+h, Yop) 
Fig. 1 


Further, in consideration of the manner in which V, and W,,; enter (cf. (4.1)), 
it is only natural to require 


v%,20 and w,=W,,;+W,2;20 HAs fate 


This simply says that when an increment K,, for U is found, U, is altered “‘in 
the direction of” K,;. Note that w, 20 does not imply that both W,,; and W,,, 
are non-negative. Indeed, in the system of parameter values given in the next 
section, Woi.4—= — 1 and Wy44=2 so that wy= 1. This is reasonable to permit, for, 
in (4.1), this means that Us$4 is given by 


Use.a = Po + Koea + (Kto4 — K214) + By. 


Thus , is altered principally by the amount K%.4, with an additional alteration 
equal to the amount the second increment K¥»4 exceeds the first increment EG an 
(This is reminiscent of the over-relaxation technique common in numerical work.) 

The determination of the parameter values of §10 was carried out under the 
further requirement that all the A,;; be non-negative. An argument such as 
that just given would imply that this requirement is perhaps too stringent. 
From the preceding extra requirements that v,, w, and A,;, be non-negative 
it follows that the intermediate parameters a; and B,; are noneiceine In 
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carrying out the actual determination of a system of parameters, it was also 
assumed that y,; and 6,; were to be non-negative. 


It is seen from the preceding sections that for each case, uw, p and q, there 
are 27 unknown parameters and 11 requirements. The additional requirements 
just discussed reduce somewhat the apparent underdeterminateness of the system. 
However, the non-linearity of the equations still makes study of the system 
difficult. 

An important question for the future improvement of the procedure is per- 
tinent here. With the parameter values given in §10, f is evaluated fifteen times 
for each subrectangle. Could values for e, g,, v,, w,, Wyo, and W,,; be found 
which reduce the number of points of evaluation of f? With regard to the 
possible relaxation of some of the extra requirements above, the manner in 
which the parameters enter in the convergence proof in §12 should be considered. 


§ 10. Values for the Parameters 

To obtain the following system of parameter values, all of the requirements 
of the preceding sections were imposed. It was found that, in order to have 
the same values of e, used in all three cases, wu, # and q, and still yield non- 
negative solutions of (R 4)—(R 10), it was necessary that e,=4 and e,=}3 (or 
vice versa). This determined separate values 2; for each of the three cases. 
The definitions of «,; and £,; (as well as of D,, and Dj») suggested choosing 
v,—=w,=1 for all three cases. These definitions, with the known values of f%;, 
and equations (R 2), (R 3) led to fixing the values of g?. From equations (R 7) 
to (R10), y?; and 6%; could then be determined. The preceding choices and 
the interrelation between the «,; and f,; already insured the satisfaction of 


(R 11). In (R 1), values 4,,, were still free to be chosen to satisfy this requirement. 
Throughout the above analysis, a preliminary study of the system 


a+btc=d,, 
atbet - ott Be 
Cc 

a+be 4+ oie de 
was of considerable aid. The solutions a, b,c of this system can be obtained 
explicitly in terms of e, d,,d, and d,. Equations (R 4), (R$) and (R6) form 
a system of this type, while (R 7), (R 8) and (R 9), (R10) form related systems 

be+-—=d,—a. 


Here, the first unknown, a, is suggested as a parameter in this system. 


Let O and I denote the zero and unit matrices, respectively. The system 
of values obtained using the analysis outlined above is the following: 


In all, three.cases,; 4,1) 2g 45. = 3 and Kah hal 
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For the computations for w: 


U ie 2 Ewe 
£1 =S5| 52a & 9 
we = Wy = diag (1, 1,1, —1, 1) a1 = diag (1, 1,1, 1, —1) 
W=1 Wy, = diag (0, 0, 0, 2, 0) W,', = diag (0, 0, 0, 0, 2) 
1 == 
Ati a5 Axis 0 A313 = 0 
Ajo3=0 Ago3=0 Ass , 
1 38 ae 
Aisa = Ag33= 28 A3s Behe sine 


For the computations for p: 


All other A;;4 are zero. 


Note that no computations with A=3 are needed in the case for # because 
Azg;4=0 (i=1, 2,3). Hence, Wi, and W are not given above. 
For the computations for g: 


ei=1 e341 8 =3 
Wii = Wei We = Wi 
Wi, = We W3h = Wes 
Ai35=f As35=4- 


All other A,;, are zero. 


Note that no computations with A= 2 are needed for case g because Ay; ;=0 


(G46 ORS y 


The values for the A,;; given above can be substituted into the formula 
(4.3) for A. This yields a convenient direct computational formula for A; no 
reference to the actual values of the A,,;; is necessary in the applications of 


R-K 2. Thus (cf. also equations (12.2)) 
h 
ph 
(10.1)  A=| ph? [sof (Uy) + sf (Urs) + Saf (Use) + Sef (Us)] | +B. 
ph lat (Uf) + (UP 
\ haf (Ute) + $f (Us) | 


To obtain this result, use has been made of equations (4.1) and (3.5) as well 
as the appropriate values of e, and g?. 
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§ 11. Determination of U, and B? 

There are a few additional considerations necessary for carrying out the 
complete numerical process on the whole rectangle OX x<a, OX y<b. The first 
of these concerns the determination of the initial vector U, for subrectangles 
along the boundaries x=0 and y=0; the others concern the determination of 
Sate t’ and t”. 

For simplicity, it will be assumed here that the mesh is given by &,=kh 
and 4,=lph where h=a/m and yh=b/n for m,n positive integers. The sub- 
rectangle for which %»=é,, yy=7, will be denoted by R,,. 

If both k>0 and/>0, then U, is given by the values of uw, and g at (kh, 1h) 
found by computations on R,_; )-1. 

Suppose k=0. (The considerations for /=0 are similar.) In this case, one 
can set u=t(lph) and q=t'(lph). To obtain f,, one can make use of the 
second of equations (1.5). Indeed, setting 


p(r,p) =1(0,7, T(r), b, v'(7)), 


y 


(11.1) p (0, y) =0'(0) + Jy(r 600, r)) dr. 


one has 


From this, f)>=£(0,/mh) can be found by the ordinary Runge-Kutta process. 
In one form of this process, it is necessary to know the values of t and t’ not 
only at the mesh points 7mh, but also at the intermediate points (+ 3) yh 
(j=0,1,...,~—1). Another form of the process uses, instead of the latter, 
the values of t and 1’ at the points (7+ $) pf and (7+ 2) ph. This form is perhaps 
somewhat preferable since the same values are used to find 7’ below. 


The quantities $’, t’, 8” and f”’ are needed in the definition (3.10) of BS. 
They are defined to begin with by equations (3.8). For the subrectangle R,, 
under consideration, recall that s’(%.) means #(&,%,), for example. In the 
numerical procedure the values of s(%)+2h), t(vp t3@A), etc., are not known; 
however, these can be approximated via Taylor expansions. Equations (3.11), 


(5.15) and (5.16) —(5.18) show that the resulting approximate values for 8’, ¢’, $”” 


and ¢”” should be accurate through order h. 


tr 


Once found, the formulas for computing S$’, t’, 5 and ¢’’ then become in 
actuality the defining formulas for these quantities; the equations (3.8) must, 
for the convergence proof of §12, be considered as merely motivation. Among 
the possible formulas, those given below have been found suitable for use in 
the convergence proof, and maintain the prescribed accuracy. As a word of 
caution, it is not clear that for other formulas the convergence proof is valid. 
For example, if §’’ is found using w instead of #, then one is led to simultaneous 
inequalities. Probably §’’ may be computed via any method of differences in #. 

In keeping with the view that the equations below constitute new definitions, 
equality signs will be used for the most part instead of ~ or an added O(h?) 
to indicate that the quantities of (3.8) are given approximately. 

From the original equations (3.8) for s’ and t’, and from a few terms in the 
Taylor series and simple forward finite difference substitutions for s’’(x)) and 
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t’’(Vo), one finds on R,; that 


§ = 28'(%) + 38 (% + 2) = $P (Ex, m) +30 (Env) » 
(14-2) F300) + 0 + A) = Balen) +39 oe tea) 
These equations apply to all subrectangles. x 
Since in (3.10) the products hs’ and wht’ are used, it is these that will be 
found here. From (3.8) 
(41.3) hg" = 3 [s'(%) + $A) — 8'(%)]. 


The formula to be given for 1” is to agree with this through order h?, Expanding 
s'(% + 2h) in a Taylor series about % and replacing s’’(%9) and s’’’(%9) by central 
finite difference formulas, neglecting terms of order h’, one obtains the result 


(11.4) hs" =2[58'(% +h) — 48’(%) — s'(xo — fy]. 
Similarly, 
(11.5) ght’ =§ [St (yo+ ph) — 4t(¥o) — to — @h)]. 


These equations can be used whenever k>O and />0, respectively. 
When k=O and /=0, if o’(2/) and t'(2yhA) are not known, then one can use 
the formulas 


(11.6) hé”" = 3, [o'(3h) — 70'(2h) + 350"(h) — 290"(0)], 
(11.7) phi’ = sq [t'(3ph) —70'(2ph) + 350'(ph) — 297'(0)]. 


The equation for 4s’’ when k=O and />0 is again derived using Taylor 
expansions. Indeed, 


(11.8) s’(2-h) = s'(0) + = hs'"(0) + : “ s""(0) + O(h'), 
and 
(11.9) s‘(3 h| Sh) ; hs!’(h) + 5 oe s'’"(h) + O(h3). 


Into (11.8) one may substitute 
s’"(0) = s"(h) — hs'"(h) + Oh?) and s’”"(0) =s'"(h) + O(h). 
Combining the result with (11.9) so as to eliminate the term in s’’"(h), one has 
s'(2h) = 4[s'(0) + 85'(h) — 2h s'"(h)] + O(R3). 


Here, s’’(h) may be replaced by #,(h, yg). Then (11.3), through terms of order 
h*, yields 

(11.10) AS" = $[s’(h) — s'(0)] — hp, (A, yo). 

In this, s’(h)=p(h, yo) and s‘(0)=£(0, yo) are known from computations on the 
subrectangle Ry ;_1. 


To determine p,(h, yp) through order h? two methods are offered. First is 
the central difference formula 


Dy Cea) ee p(2h, dibetat Vo) + O(h2) 
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which, substituted into (11.10), yields on Ry ; 
(11.11) hs" = — GP (E,m) + $b (Ex, m) — $P (0, m1). 


This uses the value #(é,,7,) from the computations on R, 7-1 as well as the 
results from Rp ;_4. 


Another method for determining #,(, Yo) is by an inductive formula: 
Px (2, 1) = be (A, M1) +392 [Dey (sm) + bey Vs M-1)] + OF). 


In this, one may use the original differential equation to replace p,, by ae A 
Thus * 
(14.12) p.(,m) =Prllb, m1) +49 [FU m) —F (0m) +F (hs 2) —F 0, 71-2)] +O), 
in which F(x, y) means f(x, y, u(x, y), £(%, y), ¢(% y)). This does not use any 
results from computations on R,;_,. The resulting formula for §’’ may there- 
fore be somewhat more desirable as regards the overall computational plan for 
R-K 2. 

In an analogous manner, in case /=0, one obtains the alternative formulas 


(11.13) ght = —44(Ex,s) + $4 (Em) — 39 (&, 0) 
and, if also k>0O, 

(11.14) pht’=$(9(E.m) — 7&0] — $9 2% (Ee. m) 
in which 


Gy (E, YP h) = Vy (Shen Pp h) a5 


(11.15) 
aia se Lee, pi) = ELE, 0) FE, Gh) — (eee, 0)]. 
§ 12. Convergence of the Numerical Approximations Obtained by R-K2 

Here we shall use the notation and terminology defined in [3] since the Theorem 
below is proved by applying I-Theorem 3%. 

The numerical process R-K 2 defines a function from the nodes of a mesh 
on R into 3-dimensional real space ¥': to every node point (&,,7)) is attached 
the triplet of values (u(&;, 7), (Ee, m), 7(&%.)). This function is clearly a vine 
as defined in I-§4. The following theorem establishes that w= w(&,,7;) is in fact 
an approximation to the solution of 2-HP. 

It is to be recalled that the regularity assumptions imposed on f in §1 are 
sufficient to insure the uniqueness of the solution of 2-HP. 


Theorem. Let /, o and 1 satisfy the regularity assumptions of §1. If a con- 
vergent sequence of meshes be considered in R, and if for each mesh a vine be found 
by means of R-K 2 using the parameter values of § 10, then the sequence of vines 
so determined converges uniformly to the solution of 2-HP. 


Notation. In the foregoing sections, computations on only one subrectangle 
R,, have been considered. In the following, further subscripts (, /) will be used 
to distinguish quantities used on each R,,, for example Uj, U;, §’". Also, u; ; 
will denote u(E,, 7); Dp. and g,,; are defined similarly. ©? 2 2 

Proof. In general terms, to prove the assertion, according to I-Theorem 3, 
two things are to be shown: (1) For each mesh the node function determined 
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by R-K 2 can be represented in the form of equations I-(4.1) (recall the remark 
preceding I-Theorem 3): 


(412.4 a) Uh +1, pal, = =Op4 Ant Ti44 — 09 A yy 2 Das "2 Ak, A NB» 
a 

(12.4 b) Phtt,141 = Ona =a x Cie ANp 

(12.4 ¢) Geeti4a = Tia t+ Bhi de. 


(2) The quantities C’, are pieane as in I-§4 by functions @ (x, y, u, , g) 
which converge uniformly to f(x, y, 4, ~,g) on Rx¥Y and which are constant 
on each cell of a convergent mesh .@*” on Rx¥Y. I-Theorem 3 can then be 
applied. In the present case, because the solution of 2-HP is unique, the whole 
sequence of vines (not just a subsequence) must converge to the solution. 

The superscript » will be suppressed throughout the following. Further, for 
simplicity, the mesh lines in R will be supposed equally spaced: £,=kh, 4,=lph 
where h=a/m and ph=b/n (m, n positive integers). Thus 4&,=h and Ay, =¢h. 
The convergence of the sequence of meshes as »->oo is assured if h->0, since 
for R-K 2 it is assumed that is bounded (see §1). 


To show (1), one finds from equations (4.4), (10.1) and (3.4) that 


(12.28) Upp rps = Meta, + Metta — Mea + 


+ ph las HU, Oy) + a 1(U 2) + oa (Us #2) + 64 (C52) 


(12.2b) Perrita = Petia + Vh 100 i hoe ; 


(k, !) 


(12.2c) ein r41=%n141 +h [a (02) HE a 1(Cs)) : 
, > ) 


The three quantities C?,, Ci,1,; and Cz,,, are then defined as the coefficients 
of ph, yh and h, respectively, in these three equations. (Note also the subscripts 
k+1 and /+1 in the right members of (12.2b) and (12.2c), respectively; they 
are due to the matrix B in (10.1).) 

It is then easy to show by a double induction on & and / (using the natural 
partial ordering of nodes, cf. I-§3) that equation (12.2a) implies that u,. 114118 
given by an equation of the form (12.1a). Similarly, by induction, p and q are 
given by equations of the form (12.1b) and (12.1c). 

From the form of equations (12.1) and the definitions of the Cas it follows 
from the boundedness of f that wu, p and g are bounded (cf. I-§4). 

The second thing which is to be shown to prove the Theorem concerns the 
convergence of the C},; to certain values of /. More specifically, Ci, ; is to converge 
uniformly to /(&, mj, Up,1, Prt» Yx,1) AS h->O (cf. the statement of Lemma 2 and 
the definition of C,, ; preceeding it in I-§4). From equations (12.2) one obtains 


the formulas 7 s oy a ae ” 
Cir= gat (Up) +47 (is) + 22 ¢ (Us) + 227 (OS), 
(k,l) (k, 1) (k, l) (k, 1) 
(12.3) Cree = i/ (U4) + 37 (U8), 
(k,l) (k,l) 
Giana a af (US) = 3 f (Us) : 
(k,l) (k,l) 
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Note particularly the subscripts k+1 and/+1 in the last two of these equations. 


Recall that U) is the vector notation for (Sh, Mir Ue,1» Pe,t» M,1)» Thus, it is to 
(k, 1) 
be shown that, as h->0, 


( 
F ( 
(12.4) Cras f(U)), 


Chaat U). 
(k, 1-+1) 


By the continuity of f, this will be true if 
Ujs > UY (AHMED, 38 


(k, L) (k, Ll) 
(12.5) Up; (Aeatle 2) 
(kt) (k-+1, 1) 
Ups, (A=1,3) 
(k, L) (k, 1+-1) 


uniformly in k and / as h-0. 


With this in mind, we write out explicitly the vectors Uj}. From § 10, 
Vo=Wi+ W2=I (w=, p, gq). Hence, from (4.1) and (3.9), 


E+ gah 
m+ agiph 
Pech) Ueneens S| trieeahs tats Pe : ane en ey 
ee eat the Tie bart ebhs” As ee hf (Uh) 
Gr,1 + & a : ni + gah /(Uo) 


The expression for U;5 is exactly the same with the sole exception that the 
argument of f (viz U,) is Ux. 
(2) (2) 
For Uf,, recall from §10 that gf=g=1. Also é,+h=&,,,. Hence, 


ae 
M+e ph 
(12.7) Up, =U, + Bi + Kh, =| “ne Fa pies pe ee 
Se ade. Pras phi Uy 


k,l) 


aiteaphl + hiU 0) 
,L) (k,1) 


In the last component of this vector, q,, age AtaL >) can be considered as an ap- 
1) 


proximation to g,,,) and, in fact, according mt (12.2), approaches g,4,) uni- 
formly since f is bounded. An expression almost identical with that just given for 


UP, holds also for Uf, the only change being that the argument of / becomes U/,. 
(k,1) 
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Similar expressions may be found for Uj;. 

In order to show (12.5), one sees from (12.6) and (12.7) that it is enough to 
show that hs’, ght, hs’ and ght" all tend toward zero with h. Indeed, since 
|{|<M, the terms in f cause no trouble. From equations (11.2) defining (for 
the purposes of the numerical computation) the quantities §” and 7’, it is clear 
that h$’—>0 and pht' >0 as h->O if s’ and ¢’ are bounded, that is if and q 
are bounded; but this has already been pointed out above. 

From equations (11.4) and (11.11) defining 4S’, if it is shown that 
\Pr41,141—Px,141| tends to zero uniformly in k and/ as h-+0, then it follows that 
h&’->0 uniformly in k and / as AO. (If the alternative formula involving 
(11.12) is used for 4S’ when k=O, difference estimates much like those below 
yield the same conclusion.) Since the considerations for yht'’ are completely 
analogous, they will not be given here. 

From (12.1 b) 


l 
(12.8) | Prta,14 = Pr, i+a| = | x41 = on a 2 | Ch41,8 > Ci, a| Ang. 


In this, by (12.3) and because the first derivatives of / are bounded, say by i, 
one has 


(12.9) | Chita — Cie] SaL| UR — — Ub] +22 | Ub — Ul. 
(k,B) (k—- (k,B)  (k—1,8) 


The vectors Uf, are given by (12.7) with the argument of f changed to U?,. From 
(US, Tighe A==41, D, 


Uf, — Uf 
(12.10) (k,B)  (k—-1,6) h 
0) 
Uni1p— Unpt a— h( i = {' ) +e, ph? [{ (UA) =F UA | 
(k,B) (k—-1,8) (Rk, B) (k—1, B) 
Prit,a — Pep + ep h(t(UA) —f( UA )] 
= Ge) (k—1, B) 
Vk, 8 — In—1,8 + rp h( (oe be ) + h(t (Uf) — f( Uf, )] 
(k,B) (Rk—-1,8B) (k, B) (k—1, B) 


From equations (12.1) and the boundedness of f, the differences 


(12.11) [41,8 — %a,p| and [de 9 —%-s,0l 
approach zero uniformly in k and 6 as h->0. Because f and t’ are bounded, 
the terms in (12.10) involving these appr oach zero as h->0. 


From equations (11.5) defining wht” when k>O and />0, 
(k,l) 


itt ah jit - 5 
a Peal = 6 | 90,84 aan Tn—1,8-41| =e ace — In—1, | +4 | Vk, B—1 ar Tk--1,p-1 ‘ 


The same estimates as those concerning (12.11) thus show that 


ph| jc ae | 
(k,B) (k-1,6) 
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approaches zero uniformly in k and B as h->+0. Similar statements apply when 
equations (11.13) are used for pht’’. 


These considerations show that, given ¢,, there is a 6, such that if h<6, 
then in (12.9) 


[Craa¢— Cie) Sa4+ L| pess.p — Pral- 


Thus, in (12.8), given e>0 there is a 6 such that <6 implies 


1 
| Pes ff ei ee Pers S64- E 2. Prti,6 — Pr, al Ang. 


Here, the fact that o’(x) is uniformly continuous has been used. By applying 
D1az’s Lemma (I-Lemma 1) it is then seen that 


lPraai41 — Paral Sle+ L| fr41,0— Prol Anoter’- 


Since p; 9>=0' (§;) and o’(x) is uniformly continuous, it follows that |P.44 144—Pp 241! 
—>0 uniformly as h->0. This establishes the desired convergence in (12.5) and 
hence that in (42.4). 

In order to define the functions @'(x, y, u, p,q), let MW be the given mesh 
on R, 4, a mesh on ¥, and.@* the product mesh. .w x4 on RXV (cf. I-§4). 
On the cell of .@* containing the point U, let 

(k,l) 


G! (x,y, u, 2, g) =Ch, (ie ep Mie eee MUS OS A yee F0) c 


On each remaining cell let a value taken on by f in that cell be chosen. Let 
the value of @ on that cell be the chosen /-value. The convergence expressed 
by (12.4) then assures that the functions @ converge uniformly to /. This com- 
pletes the proof of the Theorem. 
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The Invariants of Six Symmetric 3= 3 Matrices 
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1. Introduction 


In a previous paper [/] it was shown that every irreducible invariant, under 
the orthogonal transformation group, of any number of symmetric 3 x3 matrices 
is necessarily of degree less than or equal to six in the elements of the matrices. 
In [7] and [2] the invariants* of five or fewer symmetric 3 x3 matrices have 
been considered, and finite integrity bases have been given for these cases. 

This paper will consider the invariants of six symmetric 3 x3 matrices, 
which will be denoted a, b, c,d, e,f. As shown in [/], an integrity basis for 
these matrices consists of the integrity bases for the six sets of five matrices 
which can be selected from a, b, e, d, e, f, together with the invariant 


trabcdef (1.4) 


and invariants derived from (1.1) by permuting the matrices a, b, ec, d,e,f. Of 
the invariants which can be so formed, only a limited number are independent, 
and the purpose of this paper is to determine an integrity basis for the six 
matrices which contains the smallest possible number of elements. 

Since every irreducible invariant of any number of 3 x3 symmetric matrices 
is of degree six or less in their elements, it follows that each such invariant 
can involve at most six distinct matrices. Hence an integrity basis for any 
number of 3 x3 matrices consists of the sum of the integrity bases for the matrices 
taken six at a time in all possible combinations. Thus the results of this paper, 
with those of the previous papers mentioned above, enable integrity bases for 
any finite number of symmetric 3 x3 matrices to be constructed. It is therefore 
unnecessary to proceed to study the invariants of seven or more matrices, and 
this paper completes the analysis of the invariants of 3 x3 symmetric matrices. 


2. Relations between the invariants 


The invariants to be considered are those of the form (4.1) and forms derived 
from it by permutation of the matrices. Frequent use will be made of the follow- 
ing lemmas, which are immediate consequences of the definitions of matrix 
multiplication and the trace of a matrix. 


* Throughout this paper, invariance is understood to mean invariance with respect 


to the orthogonal group. It is immaterial whether the full or the proper orthogonal 
group is considered, 
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Lemma 1. The trace of a matrix product is unaltered by a cyclic permutation 
of the factors of the product. 


Lemma 2. The trace of a matrix product formed from symmetric matrices is 
unaltered by reversing the order of the factors in the product. 
There are sixty mvariants of six matrices a, b,c, d,e,f of the type (1.1) 


such that no two of the sixty may be equated to each other by applying Lemmas 1 
and 2. These sixty may be taken to be 


trabcdef, trabcdfe, trabcedf, trabcefd, trabcfde, trabcfed, 
trabdcef, trabdcfe, trabdecf, trabdefe, trabdfce, trabdfee, 
trabecdf, trabecfd, trabedcf, trabedfce, trabefed, trabefde, 
trabfcde, trabfced, trabfdce, trabfdec, trabfecd, trabfede, 
tracbhdef, tracbdfe, tracbedf, tracbefd, tracbfde, tracbfed, 
tracd bef, tracdbfe, tracdebf, tracdfbe, tracebdf, tracebfd, eH) 
tracedbf, tracefbd, tracfbde, tracfbed, tracfdbe, tracfebd, 
tradbecef, tradbcfe, tradbecf, tradbfce, tradcbef, tradcbfe, 
tradcebf, tradcfbe, tradebcf, tradecbf, tradfbce, tradfcbe, 
traebcdf, traebdcf, traecbdf, traecdbf, traedbecf, traedcbf. 


We use the notation 
LMxyz=xye+uzytyxetyerxt cuytzye. (2.2) 


Also, as in [J], [2] and [3], the notation J=0 will be employed to express the 
fact that an invariant J is reducible. It was shown in [3] that, if w, x, y, 2 are 
3 x3 matrices, not equal to the unit matrix J, then 


tro Lxyz=0. (2.3) 


Moreover, all the properties of invariants of 3 x3 matrices which were obtained 
in [1], [2] and [3] may be derived from Lemmas 1 and 2 and relations of the 
form (2.3). 

We will consider all possible relations between the invariants (2.1) which 
can be obtained by making substitutions in (2.3). First, it is noted that if 
w, 2, Y, z are any 3X3 matrices (not necessarily symmetric), then 


trw Lxyez=tradLyzw=trylLzwe=treLwery. (2.4) 


These relations are readily verified by introducing relations of the type (2.2) 
into (2.4) and applying Lemma 1. 

In order that (2.3) shall represent a relation between the invariants (2.1), 
way z must be a matrix product which contains each of the factors a, b, c,d, e,f 
once and once only, and each of w, x, y,2 must contain at least one of the 
factors a, b, ec, d, e, f. There are two possibilities. 

(i) One of the products w, x, y, # is a product of three of the factors a, b, 
c,d,e,f, the remaining three products from w, x, y, each consisting of a 
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single factor from a, b, ¢, d, e, f. From (2.4) it follows that it may be assumed 
without loss of generality that the product which contains three factors is the 
one which precedes the summation sign in (2.3), so that (2.3) takes the form 


tr (abe) Sdef=0, (2.5) 


or a form obtained from this by permuting the matrices a, b, ¢, d, e, f. 

(ii) Two of the products w, «, y, # are products of two factors chosen from 
a, b, c,d, e, f, the remaining two products each consisting of a single one of 
these factors. From (2.4), it may be assumed without loss of generality that 
the product preceding the summation sign in (2.3) is one of the products which 
contains two factors, so that (2.3) takes the form 


tr (ab) XY (ed) ef =0, (2.6) 


or a form obtained from this by permuting the matrices a, D, ¢, d, e, f. 

The relations (2.5) and (2.6), and relations derived from them by permuting 
the matrices, include all possible relations of the type (2.3) between the invariants 
(2.1). Since all the relations between invariants which can be obtained from 
results given in [1], [2] and [3] may ultimately be based on (2.3) and Lemmas 1 
and 2, it is sufficient to consider relations of the types (2.5) and (2.6), in con- 
junction with Lemmas 1 and 2, in order to determine how many of the invariants 
(2.1) are independent. It is, however, convenient to use some further relations 
which, from the manner in which they are derived, must be expressible as linear 
combinations of relations of the forms (2.5) and (2.6), although we shall not so 
express them explicitly. In [3] it was shown that, if w, a, y, 2 are defined as 
in (2.3), then 

traye2w =0, (2.7) 
and that if x, y, 2 are also symmetric, 
trae? y?2? = 0. (2.8) 


In (2.7), replace x, y, , w by a+b, e, c+d, f respectively. Then, on expanding 
and using the relations (of type (2.7)) 


trWece’f=0, tr@e@f=0, trbec*f=0, trb¢ed’f=0, 


equation (2.7) becomes 
tr(ab + ba) e(cd+de)f=0. (2.9) 


In a somewhat similar way, we obtain from (2.8) the relation 
tr (ab + ba) (ed + de) (ef + fe) =0. (2.10) 

We will use the notation 
A(a, b,c; d,e,f)=trabe2 def, (2.14) 
Bia, b; ¢,d:e,f) =trabZ (cd) ef, (2AD) 
C (a, b; €,d:e, f) r(ab+ ba) e(cd+de)f, (2.13) 
) (2.44) 


=t 
D(a, b; ¢,d; e,f) = tr(ab + ba) (ed + de) (ef+fe). 
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From the definitions (2.11), (2.12), (2.43) and (2.14), and Lemmas 1 and 2, the 
following relations follow immediately. 


(i) A(a, b,c; d,e,f) =A(e,b,a;d,e,f), (2:45) 
and A(a, b,c; d,e,f) is unaltered by permuting d, e, f in any manner. 


(11) B(a,b; ¢,d:e,f) = B(e,d; a, b:e, f) 


(2.16) 
= B(b,a;d,c:e,f) =B(d,c; b,a:e,f), 
and B(a, b; c,d: e,f) is unaltered by interchanging e and f. 
(iii) C(a, 65 c,d: ¢€,f) = Cle, d; a, 0: e,f), (247) 


and C(a, b; e, d: e, f) is unaltered by interchanging a and b, or e and d, or e and f. 

(iv) D(a, b; c,d; e,f) = D(a, 6; €, f;¢,d) =D(e,d; a, b; e,f) : 

Dl 
= D(e,d;e,f;a,b)=D(e,f; a,b; c,d) =Dle,f;¢,d;a, b),| ) 

and D(a, b; ¢, d; e, f) is unaltered by interchanging a and b, or ce and d, or e and f. 

Taking into account these symmetry properties of A(a, b,c; d,e, f), 
Bia, b; e,d:e,f), Cla, b;e¢,d:e,f) and D(a, b; e,d;e,f), there remain for 
consideration 60 relations of the type A(a, b, c; d, e, f) =0, 90 relations of the 
type B(a, b; c,d: e,f) =0, 45 relations of the type C(a, b; ¢,d:e,f)=0, and 
15 relations of the type D(a, b; c,d; e,f)=0. These relations are, of course, 
not all independent. 

3. Expressions of the type A(a, b, c; d, e,/f) 

In this section it will be shown that the expression A(a, b,c; d,e,f) can 
be expressed as a linear combination of expressions of the type C(a, b; ¢,d: e, f). 
We have, from (2.11) and (2.15), and Lemma 2, 

2A(a, b,c; d,e,f) =tr(abe 4+ cba) Ldef 

= tr{a(be + cb) — b(ac + ca) + c(ab + ba)} def 
= tr{a(be + cb) — b(ac + ca) + c(ab + ba)}x 
x {d(ef+fe)+e(df+fd) +f(de+ed)}, 

from (2.2). Hence, by (2.13) and Lemma 1 


2A(a, b,c; d,e,f)=C(b, c; e,f:a,d)+C(b,c; d,f:a,e)+C(b,c; d,e:a,f)— 
— C(a,c; e,f: bd) — C(a,c; d,f:b,e) — C(a,e;d,e:6,f)+ (3.4) 
+ C(a, b; e,f:c¢,d)4+ C(a,b; d,f:c,e)+C(a,b;d,e:e,f). 


Thus the relations of the form C(a,b;c,d:e,f)=0 imply the relations 
A(a, b,c; d, e,f)=0, and the latter relations need not be considered further. 


4, Expressions of the type B(a, b; c,d: e,f) 
We now consider expressions of the form B(a, b;c,d:e,f). It may be 
verified by using Lemmas 1 and 2, with (2.2), (2.11), (2.12), (2.13) and (2.14), 
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that B(a,b;¢,d:e,f)+ Bac; b,d:e,f) 
— A(d,f,a;b,c,e)+A(d,e,a; b,¢,f) = 

<= A(b, ¢, a; d;e,f) —Alb, Gd; a, ef) = 

— A(b,c,e;a,d,f) — A(b,¢,f; 4,d,e)+ 
+ 2C(b,¢;e,f:a,d)+ C(a,d; e,f:b, c) + C(a,d; b,c: e,f) + 
Fils D(a, d; b, Cc: e,f). 
By relations of the form (3.1), the expression on the right-hand side of (4.1) 
can be expressed as a linear combination of expressions of the types C(a, b;¢,d:e, f) 
and D (a,b; ¢,d;e,f). 

Also, it may be verified that 
B(a,b;¢,d:e,f)+ B(a, 6; d,c:e,f) 
= C(a, b; c,d: e,f) + D(a, b; ¢,d; e,f). 

From (4.1), (4.2), relations obtained by permuting the matrices a, b,e and d 
n (4.1) and (4.2), and (2.16), it follows that all of the expressions which can be 
formed from B(a,b; e¢,d:e,f) by permuting the matrices a, b, c,d can be ex- 
pressed as linear combinations of B(a, b; ¢,d:e,f) itself with expressions of 
the types C(a, b; c,d: e,f) and D(a, b; c,d; e,f). Thus of the 90 expressions 
of the form B(a, b; c,d: e,f) described in Section 2, 75 may be expressed in 
terms of the remaining 15 and expressions of the types C(a, b; c,d: e,f) and 
D(a, b; e, d; e, f), and need not be considered further. The fifteen expressions 
which are retained correspond to the fifteen ways in which two matrices can 
be chosen, from, the six matrices a,b,c, d,e,f. Thus Bia, b; c,d: 6,7) cor 
responds to the pair (e, f), and so on. The fifteen expressions which are retained 
are taken to be 


(4.1) 


(4.2) 


B(a,b;e,d:e,f), Bla,b;ef:e,d), Bi(b,e;d,e:af), 
B(a,b;ec,e:d,f), B(a,e;d,e:b,fy, B(b,ce;d,f:a,e), 
Bia, b;¢,f:d,e), Baye; df psey) ) PB10; 6: e. 7: ad), (4.3) 
Bia,b; d,enc,f);. Bia,e: €,f2-b,d).\ Boyd; €,f 5 a,c), 
B(a,b; d,f:e,e), Bilagd; e,93.0; Coy |B (6, d; e, f a,b). 


It may further be verified, by using relations of the types (2.11) and (2.12), 
that B(a,b;¢,d:e,f) + B(a,b; ¢,e:d,f)+ B(a,b;c¢,f:d, e) 
= A(d,a,b;c¢,e,f) +A(e,a,b;¢,d,f)+A(f,a,b;e, e,f). 
Hence, from (4.4) and relations of the type (3.4), 


(4.4) 


Bla, b,¢,d:e,f)+B(a,b;¢,e:d,f) +B(a,b; ¢,f:d,e) (4.5) 


can be expressed as a linear combination of expressions of the type C(a,b;e,d:e,f). 


Noessentially new relations between the expressions of the form B(a,b;e,d:e, f) 
can be obtained by permuting a, b, ce in any manner in (4.5). For example, 
interchanging b and ¢ in (4.5), we find that the expression 


B(a,c;b,d:e,f) + B(a,c; b,e:d,f) + Bla,ec; b,f:d,e) (4.6) 
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can be expressed asa linear combination of expressions of the type C(a,b; ¢,d:e,f). 
This, however, merely expresses a fact that can be readily deduced from (4.5), 
expressions of the types (4.1) and (3.1), and the relation (easily verified from 
(2.11) and (2.14)) 


D(a,d; b,c; e,f)+D(a,e; b,c; d,f)+D (af; b,c; d,e) 


AN 
=A(b,c,a;d,e,f)+A(b, c,d; a,e,f)+A(b,¢,e;a,d,f)+A(b,¢, f;a,d, e) 


Again, interchanging @ and b in (4.5), it follows that the expression 


B(b,a;c,d:e,f)+ B(b,a;c,e:d,f)+8(b,a; ¢,f:d,e) (4.8) 
can be expressed as a linear combination of expressions of the type C(a,b; ¢,d:e,f). 
This expresses a fact which can be deduced from (4.5), relations of the types 
(2.16), (3.1) and (4.2), and a relation similar to (4.7). Also, using (2.16), it is 
evident that (4.4) is unaltered by permuting d, e and f in any manner. It 
follows that of the relations of the form (4.4), it is necessary to consider at most 
20 relations, these 20 corresponding to the twenty ways in which three matrices 
can be chosen from a, b, c,d, e,f. Thus (4.4) corresponds to the trio (a, b, e), 
and so on. All other relations of the form (4.4) can be derived by means of 
relations already considered from 20 relations selected in this way. The twenty 
relations retained for further consideration are taken to be 


Bia, b;c,d:e,f) + B(a,b;¢,e:d,f)+B(a,b;¢,f:d,e)=q,, 
Bia, b; d,c:e,f)+ B(a,b;d,e:e¢,f)+5(a,b; d,f:c¢, e) = gz, 
B(a, b; e,e:d,f) + Bia, b; e,d:c,f) +_B(a,b; e,f:¢,d) = 9,, 
Bia, b; f,e:d, Rs B(a,b; f,d:c,e)4+ Bia, b; f,e: c,d) =qy,, 
B(a,c;d,b:e,f) + B(a,c;d,e:b,f) + B(a,c;d,f:b,e)=¢;, 
B(a,c;e,b:d,f)+ B(a,c;e,d:6,f)+B(a,e;e,f:b,d) =, 
B(a,c;f,6:d,e)+ B(a,c;f,d:b,e)+Bia,c;f,e:b,d)=¢q,, 
Ba, eye B(a,d;e,c:b,f)+ Bia, d;e,f:b,c) =9,, 
B(a,d;f,b:e,e)+ B(a,d; f,c:b,e)+ B(a,d; f,e: b,c) =Qp, 
B(a,e;f,b:¢,d)+B(a,e;f,c:b,d)+ B(a,e; f,d: b,c) = qo, 49) 
B(b,c; d,a:e,f) + B(b,c;d,e:a,f)+B(b,¢;d,f:a,e) =, 
B(b,c;e,a:d,f)+ B(b,c;e,d:a,f)+B(b,c;,e,f:a,d)=Qp, 
B(b,¢; f,a:d,e)+ B(b,c;f,d:a,e)+B(b,c¢; f,e:a,d) =, 
B(b,d; e,a:c,f) + B(b,d;e,c:a,f) + B(b,d; e,f: a,c) =P, 
B(b,d;f,a:c,e)+ B(b,d;f,c:a,e)+B(b,d;f,e:a,c) =, 
B(b,e;f,a:¢,d)+B(b,e;f,c:a,d)+B(b,e; f,d: a,c) = ep, 
B(e,d;e,a:b,f)+B(e,d;e,b:a,f)+ B(e,d;e,f:a,b) =o, 
B(e,d; f,a:b,e)+ B(e,d;f,b:a,e)+ B(e,d; f,e:a, pes Cis: 
B(c,e;f,a:b,d)+B(c,e;f,b:a,d)4+ Bi(c,e;f,d:a,b) =p, 
B(d,e;f,a:b,c)+B(d,e;f,b:a,c)+B(d,e;f,c: 4a, b) = ga, 


where the q; are linear combinations of expressions of the type C(a, b; c,d: e, f). 
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We next express the relations (4.9) as relations between the fifteen expressions 
listed in (4.3). Using relations of the types (4.1) and (4.2), equations (4.9) may 
be written 


B(a, b;¢,d:e, )+ B(a,b;c,e:d,f)+B(a,b; ¢,f:d,e)=y, 
— B(a,b;¢,d:e,f)+ B(a,b;d,e:¢,f)+B a,.D;,d;f 1 Cx) =e, 
— B(a,b;c,e:d,f)— B(a,b; d,e:¢,f) + B(a,b;e,f:¢,d) =, 
— B(a,b; ¢,f:d,e)— B(a,b;d,f:c,e)—B(a,b,e,f:c,d)=y%, 

B(a,b;¢,d:e, )+ B(a,c;d,e:b,f)+B(a,c;d,f:b,e)=y,, 

B(a,b; c,e:d,f)— B(a,e;d,e:b,f) + B(a,c; e,f:b, d) =, 

B(a,b;¢,f:d,e)— B(a,c; d,f:b,e)— B(a,c;e,f:b,d)=y,, 

B(a,b; d,e:c,f)+ B(a,c; d,e:b,f)+ Bla,d;e, DC) =e. 

B(a,b;d,f:c,e)+ B(a,c; d,f:b,e)— B(a,d; e,f: b,c) =H», 

B(a,b;e,f:c,d)+B(a,c;e,f:b,d)+ B(a,d; e,f: b, ¢) = yy ae 
— B(a,b;e,d:e,f)+ B(b,c; d,e:a,f) + B(b,c;d,f:a,e)=yy, 
— B(a,b; c,e:d,f)— B(b,c; d,e:a,f)+ B(b,c; e,f:a,d) =yYp, 
— B(a,b; c,f:d,e)— B(b, ce; d,f: a, e)— B(b, €; €,f:a,d) =p, 
— B(a,b; d,e:c,f) + B(b,¢; d,e:a,f) + B(b,d; e,f:a,c) =y4, 
—B(a,b;d,f:c,e)+ B(b,¢; d,f:a,e)— B(b,d;e,f:a,¢) =, 
— B(a,b; e,f:¢,d)+ B(b, €; e,f: a, d) + B(b, d; e,f: a, c) = y,, 
— B(a,ec; d,e:b,f)— B(b,¢; d,e:a,f)+ Ble, d; e,f: a,b) —y,., 
— Bia, ec; d,f:b, e) — B(b, ¢; d,f:a,e) — Ble, d; e, f: a,b) =y,,, 
— B(a,c; e,f:b,d) — B(b,c; e,f:a,d) + B(e,d; e,f:a, b) = yp, 
— B(a,d; e,f:b,c)— B(b,d; e,f:a,c)— Bi(e,d; e,f: a,b) = yo, 


where the y; are linear combinations of expressions of the types C(a, b; ¢, d: e, f) 
and D(a, b; c,d: e,f). The matrix of the coefficients of the 15 expressions (4.3) 
in the 20 equations (4.10) is 


= 13a eee : =) 
= ee , 1 
1 Sly a 4 j =A 
1 —1 ee ae —1 
4. x4 4 1 4 
1 ABE TS A Re SR eee eee 
Ce oo ee, Come Sue Sem pone mae Ak (4.11) 
Pk SO Re, eae Ot eee 
ted ae: 5 4 
Sd na eed . =i 
1 ie ee 1 
29 yt hae ae Cae ‘ees 1 oxy 
sb ke Gl Ve Agon be Gl Bo gy a es Sk] ee meme 
Dian We ees. A, aa neg ta ata te A, RE = 


It is easily verified that each of the last five rows of this matrix is equal 
to a linear combination of four out of the first ten rows; for example, the eleventh 
row is obtained by adding the second row to the seventh row and then sub- 
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tracting the first and fourth rows. It is also evident from the diagonal sub- 
matrix which is formed by the elements common to the first ten rows and last 
ten columns of the matrix that the first ten rows of the matrix are linearly 
independent. Hence the matrix has rank ten, and only ten of the equations 
(4.10), which may be taken to be the last ten, are independent. These equations 
may be used to express the ten expressions in the first two columns of (4.3) 
as linear combinations of the five expressions in the third column of (4.3) and 
expressions of the types C(a,b; c,d: e,f) and D(a,b;¢,d:e,f). Thus each 
of the original 90 expressions of the form B(a, b; ¢,d:e,f) can be expressed 
as a linear combination of one or more of the expressions 


B(b,c;d,e:a,f), B(b,c;d,f:a,e), Bi(b,c;e,f:a,d), 


(4.12) 
B(b,d;e,f:a,e), Bi(e,d;e,f:a,b), 


with expressions of the types C(a, b; c¢,d:e,f) and D(a, b; ¢,d;e,f). Hence 
of the original 90 relations of the type B(a, b; c,d: e, f) =0, it is necessary 
to consider only the five relations obtained by setting each of the expressions 
(4.12) equivalent to zero, together with relations of the types C(a, b; c,d: e, f) =0 
and D(a, b; ¢, d; e, f) =0. 


5. Expressions of the form C(a, b; c,d: e, f) 
It was shown in Section 2 that it is necessary to give further consideration 
to 45 expressions of the form C(a, b; c,d: e,f), which correspond to the 45 
ways in which two pairs of matrices can be chosen from the six matrices a, b, 
c,d,e,f, the ordering of the pairs being immaterial. The expressions to be 
considered are 


C(e,d; e,f:a,b),. Cle,e;d,f-.a,b),... Cie, f;d,e: a,b), 
C(bi de, fra, ce), Clove; dj fa;e), (Cb, [-dre* a,c), 
C(b,c;e,f:a,d), C(b,e;c,f:a,d), C(b,f;¢,e:a,d), 
Obeid, ve, G0, 00-50, e).4) C(O. fe, 02 a7e); 
C(b,c;d,e:a,f), C(b,d;c¢,e:a,f), C(b,e;c¢,d:af), 
C(a,d;e,f:b,c), C(a,e;d,f:b,c), Caf;de:b,e), 
C(a,¢;e,f:0,d), - C(a,e;¢,f:b,d),  Cla,f;e¢,e:b,d), 
C(a,c;d,f:b,e), C(a,d;e¢,f:b,e), C(a, f; ¢, d: bye), (5.1) 
C(a,ce;d,e:b,f), C(a,d;c,e:b,f), C(a,e;¢,d:b,f), 


CS 


( 

C( ) 

C(a,b;d,f:c,e), C(a,d;b,f:c,e), ( 

C(a,b;d,e:c,f), C(a,d;b,e:c,f), ( 
(a,b;c¢,f:d,e), C(a,c;b,f:d,e), Caf;b,c:d,e), 

Clapb; 6, e:d,f), \G(a,e; b,e:d,f), » Cia ,e;b,c¢:d,f), 
( ) C(a,c;b,d:e,f),  C(a,d;b,c:e,f). 


a,b;e,f:e,d),  C(a,e;b,f:c,d), (a, f; b,e:e,d), 
( 


C(a,f; b,d:e,e), 


Da 


) ) 
) ) 
) a,e;b,d:e,f), 
) ) 
) ) 
) ) 
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We denote by S the sum of the sixty invariants listed in (2.1). Then, using 
(2.13) and Lemmas 1 and 2, it may be verified that 


S=C(e,d;e,f:a,b) +C(c,e;d,f:a,b)+C(e,f; d,e:a, b) + 
+ C(b, d; e,f:a,c)+C(b,e; d,f: a, c)+C(b,f; d,e:a,e) + 
L C(b, c; e, f: a, d) +C(b, €; c, f:a,d) + C(b, f; ¢,e:a,d) + 
+ C(b,e; d,f:a,e)+C(b,d; ¢,f:a,e) + C(b,f; ¢,d:a,e) + 
| C(b,c; d,e: a,f) +C(b,d; ¢,e:4,f) + C(b,e;¢,d:a,f). 


Similarly, a further five expressions for S may be obtained by cyclic permu- 
tation of the matrices a, b, ¢, d, e, f in (5.2). It is convenient to use the notation 


H(c,d, e,f; a, b)=C(e, d; e,f:a,b)+C(e, e€; d,f-a, b)+C(e,f;d,e:a,b). (5.3) 
Then (5.2) and the five similar expressions for S may be written 
S=H(c,d,e,f;a,b)+H(0,d,e,f; a,c) + H(b, c, e, f; a, d) + 

+ H(b,c,d,f; a, e) + H(b,c,d,e; a, f) 

= H(e, d, e, f; a, b) + H(a, d, e,f; b, ec) + H(a,c, e, f; 6, d) + 
+ H(a, c,d, f; b, e) + H(a, c,d, e; b,f) 
= H(b,d,e,f;a,c)+H(a,d,e,f; b,c) + H(a, b, e,f; ¢,d) 4 
+ H(a, b,d,f; ¢, e) + H(a, b, d,e; ¢, f) 

= H(b,c,e,f;a,d)+H(a,c,e,f; b,d) + H(a, b, e,f; c,d) + 
+ H(a, b, ec, f; d, e) + H(a, b,c, e; d, f) 
= H(b,c,d,f; a,e) + H(a, c,d, f; b, e) + H(a, b,d, f; ¢, e) + 
+ H(a, b, ec, f; d, e) + H(a, b,c, d; e,f) 
— H(b,c,d,e;a,f) +H(a,c¢,d,e; b,f) + H(a, b, d,e; ¢,f) + 
+ H(a, b,c, e; a, f) + H(a, b, ¢, d; e,f). 


We now eliminate S from (5.4) and obtain the following five equations between 
the expressions of the form H(e, d, e, f; a, b): 


H(b,d,e,f;,a,c)+H(b,c,e,f;a,d)+H(b,¢,d,f;a,e)+H(b,c,d,e;a,f) — 
—H(a,d,e,f;b,c)—H(a,c, e,f;b,d)—H(a,c,d,f;b, e)—H(a,e,d,e;b,f) =o, 
H(c,d,e,f;a,b)+H(b,c¢,e,f;a,d) +H(b,c,d,f;a,e)+H(b,c,d,e;a,f) — 
—H(a,d,e,f;b,c)—H(a, b,e,f;¢,d)—H(a, b,d, f;c, e)—H(a, b, d,e;e,f) =0, 
H(c,d,e,f,a,b)+H(b,d,e,f;a,c)+H(b,e,d,f;a, e)+H(b,c¢,d,e;a,f) — 
—H@ ees: b,d)—H(a, b, e, f;¢,d)—H(a, b, c,f; d, e)— HA (a, b, c,e;d, f) =0, 
H(c,d,e,f;a,b)-+-H(b,d,e,f;a,c)+H(b,¢,e,f;a,d) +H(b,c,d,e; a,f) — 
—H(a, ¢,d,f;b,e)—H(a, b,d,f; ¢,e)—H(a, b, e, f;d, e)—H(a, b,c, d;e,f)=0, 
H(c,d,e,f; a,b) +H(b,d,e,f; a,c) +H(b,¢,e,f;a,d) +H(b,¢,d,f;a, e) — 
—H(a,c,d,e; b,f)—H(a, b,d,e; ¢,f)—H(a, b, c,e;d,f)—H(a, b,c, d;e,f) =0. 


if 


(5.4) 


(5.5) 
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The matrix of the coefficients of the fifteen expressions H(e, d, e, f; a, b), 
H(b, d, e, f; a,c), ..., H(a, b, c,d; e,f) in the five equations (5.5) is 


OO Vid 0M hs (10) ee 21% OF Modamoh 049Gb 50 

(ee Ome ae terete OM Oe a ES 0.8 0 

Melee ieee 6 Om aed PO 102 foe OO 1 (5.6) 
foe elo ae oO BO eat Oe Ole ad tO. A Ola 

Jaton! peas 0. wi Cami Oneea=d 20, Os hierdie 0 od =tose te 


It may be verified that the determinant formed from the first five columns of 
this matrix is non-zero, and so the matrix has rank five. It follows that equations 
(5.5) may be used to express five of the expressions of the form H(e, d, e, f; a, b), 
which may be taken to be H(e, d, e, f; a, b), H(b, d, e, f; a, c), H(b, ¢, e, f; a, d), 
H(b, c,d, f; a, e) and H(b, ¢, d, e; a, f), in terms of the remaining ten expressions 
Ha, d, e,f; b, c), H(a, ¢, e, f; b, d),..., H(a, b, c, d; e, f). Since each of the ex- 
pressions H (e, d, e,f; a,b) is a sum of three expressions of the type C(e,d; e, f:a, b), 
and each expression of the form C(e, d; e, f: a, b) occurs in one and only one 
of the expressions H(e, d, e, f; a, b), equations (5.3) and (5.5) may be used to 
express five of the expressions C(e, d; e, f: a, b) listed in (5.1) in terms of the 
remaining forty expressions in (5.1). The five expressions which can be expressed 
in terms of the other forty can be chosen to be 


C(e,e;d f- a, b) we Gb, d: ef 7d, €) Cb, e:¢, [-a, 4); 


C(b;f; c,d =a, e),,, 1 C(b, ec; d, e:.a, J). 6-7) 
It follows that the relations which can be obtained by setting the expressions 
(5.7) equivalent to zero may be derived from the relations obtained by setting 
equivalent to zero the forty expressions which are included in (5.1) but not in 
(5.7). Hence, of the relations of the form C(e, d; e, f: a, b) =0, it is necessary 
to consider only the forty relations which are obtained by setting equivalent 
to zero the forty expressions of the type C(e, d; e, f: a, b) which are included 
in the set (5.1) but not in the set (5.7). 


6. Expressions of the type D(a, b; c,d; e, f) 
In this section we consider the fifteen expressions of the form D(a,b;¢,d;e, f), 
which correspond to the fifteen ways in which three pairs of matrices can be 
selected from the six matrices a, b, ec, d, e, f. These fifteen expressions are 


Diasb 6, de fj “Dia, beer es7d,f) 2 Da, by Cf; d, e); 
Dia c.b,d- ef),  Da,e7b,e;,d,f), Di(a,c;b,f;d,e), 
Di(a,d; b,c; e,f),  Dia,d;b,e;c,f), D(a,d;b,f;¢,e), (6.1) 
D(a,e;b,c;d,f), D(a,e;b,d;c,f), D(a,e;6,f; c,d), 
Dia f;iby ce; dze)s! 1D (a,f70, dsc, 6), © 'D (asf; e; c,d). 
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It may be verified, by means of relations of types (2.41) and (2. 
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D(a,b; ¢,d;e,f) +D(a,b;¢,e;d,f)+D(a,b; ef; d, e) 


=A (a,b,c; d,e,f)+A(a,b,d;¢,e,f) +A (a, b,e;c,d,f)+ 


941s; 


14), that 


A(a,b,f; e,d,e) (6.2) 


where, by means of relations of the type (3.1), y%; can be expressed as a com- 
bination of expressions of the form C(a, b; c,d: e,f). Fifteen relations of the 
type (6.2) may be formulated, corresponding to the fifteen ways in which a 
pair of matrices can be chosen from the six matrices a, b, c, d, e, f. The relations 


is D(a,b; ¢,d; e,f) + D(a,b; ¢,e;d,f)+D 
D(a, ce; b, d; e,f) + D(a, €; b, e; d,f) + 
D(a, d; b,c; e,f) + D(a, d; b,e; ¢,f) + 
D(a,e; b,c; d,f) + D(a,e; b,d; ¢,f) + 
D(a,f; b,c; d,e)+D(a,f; b,d;c,e) + 
D(a,d; b,c; e,f) + D(a, e; b,c; d,f) + 
D(a, ec; b,d; e,f) + D(a,e; b,d; ¢,f) + 
D(a,c; b,e; d,f) +D(a,d; b,e; ¢,f) + 
D(a, ec; b,f; d,e) + D(a, d; b,f; ¢, e) + 
D(a, b; ¢,d;e,f) +D(a,e; b,f; ¢,d) + 
D(a, b; c¢,e; d,f) + D(a,d; b,f; ¢,e) + 
D(a, b; ¢,f; d, e) + D(a, d; b,e; ¢,f) + 
D(a, b; ¢, f; d,e) + D(a, ec; b,f; d,e) + 
D(a, b; c,e; d,f) + D(a, c; b,e; d,f) + 
D(a, b; c,d; e,f) +D(a,e; b,d; e,f) + 


where each of the y; is a sum of expressions of the form C(a, b; ¢, d: e, f). 


(anos oj fpdsey = 75 
Dia, ec; 6, f; d, e) = 75. 
D(a, d; b,J; ¢, e) = 73; 
Diase? Gf? 0) d) 3744; 
Di{a,f: b, €-e,d)\—=¥,, 
D (a; f; boc; de) = 75s 
Dia; 0a; ¢,e—~,. 
D (a; f;-b; €;, ©; d) =4,, 
D (a, e; b, f: c,d) =y,. 
D (a, f; b, €; ¢; d) = 749; 
D(a; f;b;, d; c, €)=y7,,, 

D (a, e;,b,.; €,f) = %4 
D (a, f; 0, e; d, e) =; 

D(a,e; b, ce; d,f) =a, 

D(a,d; b,€;e,f) = ys. 


(6.3) 


We now denote the expression on the left-hand side of the first of equations 
(6.3) by X,, the expression on the left-hand side of the second of (6.3) by X5, 


and so on to X4;. 


ox 


Then it may be verified that 


= — Xe— X,— Xg— Xg 4+ Xt Xn 4+ Xe 
= —X_+ X,+ Xg4+ Xy— Xyy— Xy 

te Xe Xo Xe Xo Rag h a 

= XetX,—X,+X,4+ Xy— XyA4 

= Xe+ X_+ Xg—Xo+ Xo + Xy— 


X43+ = (an X45, 


| 
Xie I X43 Xi4 X45; 


t Xya— Xgg— Aya t Xs, 


Xp X43 Xi X45, 
Xye+ X43— Mag X45. 


Hence, regarding (6.3) as a set of equations for the expressions D(a, b; c,d: e,f), 
the first five of equations (6.3) are dependent upon the remaining ten, and need 
not be considered further. The matrix of the coefficients of the fifteen expressions 
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listed in (6.1) in the last ten of equations (6.3) is 


jo) 


3 
6) 
O 
1 
0) 
O 
6) 
0 
1 
O 
O 
1 
O 
O 
O 
O 


ee SS SS ea OS Se eS) = 
eS Se eS SS Se =| SoS Se © 
rs GaSe Oe ee SS Soa © Se, 2] 

Se (eS eo Sea) a SS] SSS 
— 1 ey Te SS) OO) Se OS) Se 
eS. Seer 2) Se So 2 SoS) "SS 
Se CS SS See a Se iS Ss 
Soe ere Ss SS] Qe are eae 
ee oo ea 2 eyo = So S22 = 


It can be verified that the determinant of the matrix formed by taking rows 
1, 3,5, 6, 7,9, 10, 11,14 and 15 of the matrix (6.5) is non-zero. Hence the 
matrix (6.5) has rank ten, and it follows that equations (6.3) may be used to 
express ten of the expressions (6.1) as linear combinations of the remaining 
five of these expressions and expressions of the form C(a, b;¢,d:e, f). 
particular, (6.3) may be used to express 

D(a,b;e¢,d;e,f), D(a,b;e¢,f;d,e), Di(a,e;b,e;df), 

D(a,c;b,f;d,e), Di(a,d;b,c;e,f), Di(a,d;b,f;¢,e), 

D(a,e;6,e;d,f), Di(a,e;b;d;¢,f), Di(a,f;o,d;¢,e), 

D(a,f; b,e; ¢,d), 


(6.6) 


as linear combinations of 
D(a,b;¢,e;d,f), D(a,ce;b,d;ef), Di(a,d;b,e; cf), 


(6.7) 
D(a,e;b,f;e¢,d), Di(a,f;b,c;d,e), 


and expressions of the form C(a, b; e,d:e,f). Thus, in place of the fifteen 
relations of the type D(a, b; ¢,d; e, f) =0, it is sufficient to consider only the 
five relations obtained by setting the five expressions (6.7) equivalent to zero, 
together with relations of the form C(a, b; ¢, d: e, f) =0. 


7. A set of independent invariants 
It has been shown in Sections 3—6 that of the relations of the forms 


A(a, b, ¢; d,e,f) =0, (721) 
B(a,b; ¢,d:e,f) =0, (7.2) 
C (a,b; ¢,d:e,f)=0, (7.3) 
D(a,b; ¢,d;e,f) =0, (7.4) 
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it is sufficient to retain 

(i) Five relations of the form (7.2), which may be taken to be the five relations 
obtained by setting equivalent to zero the five expressions (4.12) ; 

(ii) Forty relations of the form (7.3), which may be taken to be the forty 
relations obtained by setting equivalent to zero the forty expressions which 
appear in (5.1) but not in (5.7); 

(iii) Five relations of the type (7.4), which may be taken to be the five 
relations obtained by setting equivalent to zero five expressions (6.7). 

All other relations of the types (7.1) to (7.4) may be regarded as consequences of 
the relations defined in (i), (ii) and (iii) above. Moreover, (7.1) to (7.4) include all 
possible relations between the invariants (2.1) which can be derived from (2.3): 

It remains to be seen how many of the fifty relations defined above in (i), 
(ii) and (iii) are independent. In order to do this we form the matrix of the 
coefficients of the sixty invariants (2.1) in the fifty relations. This matrix is 
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In this matrix, the first forty rows are formed from the coefficients of the sixty 
invariants (2.1) in the forty relations defined in (ii) above; the next five rows 
are formed from the coefficients of the invariants (2.1) in the five relations 
defined in (i) above; and the last five rows are formed from the coefficients of 
the invariants (2.1) in the five relations defined in (iii) above. 


It has been verified* that the matrix (7.5) has rank fifty, and that, for 
example, the determinant of the 50 x 50 submatrix obtained by omitting columns 


AD eae 50, 953,. 54, 50857, 55, 59, 00; (7.6) 


from the matrix (7.5) is non-zero. Hence the 50 relations defined in (i), (ii) 
and (iii) above are all independent, and may be used to express fifty of the 
invariants (2.1) in terms of the remaining ten invariants (2.1), and only these 
remaining ten invariants need be retained in the integrity basis. The ten 
invariants to be retained may be chosen in many ways; for example, they may 
be taken to be the invariants which correspond to the columns of (7.5) which 
are listed in (7.6), these invariants being 


tracfebd, tradcbfe, tradcfbe, tradfbce,  tradfebe, 


traebdcf, traecbdf, traecdbf, traedbecf, traedcbf. aD) 
An integrity basis for the six 3 x3 symmetric matrices a, b, c, d, e, f therefore 
consists of the integrity bases for the matrices taken five at a time, together 
with the invariants (7.7). Since all possible relations of the type (2.3) between 
the six matrices have been considered, the number of invariants from (2.1) to 
be retained in the integrity basis cannot be further reduced by the use of relations 
of this type, or relations which can be derived from relations of this type. No 
other polynomial relations between the invariants (2.1) are known to exist. 
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Bemerkung xum geometrischen Grundgesetz 
der allgemeinen Kontinuumstheorie der Versetzungen 


und Eigenspannungen 


EKKEHART KRONER 


Vorgelegt von J. MEIXNER 


In der vorliegenden Bemerkung, die als Erganzung zu einer kiirzlich erschie- 
nenen Arbeit [7] aufzufassen ist, soll begriindet werden, daB die Einsteinschen 
Gleichungen J™/= B’’ eine vollstandige Formulierung des geometrischen Grund- 
gesetzes der obengenannten Theorie darstellen. Dieses Gesetz beinhaltet die 
Forderung nach der Kontinuitat des verformten Kérpers. 

In der auf 9 funktionale Freiheitsgrade beschrankten Theorie, mit der reine 
Kontinua, d.h. solche ohne Fremd- oder Extramaterie behandelt werden, war 
das Grundgesetz frither in der linearisierten Form [2, 3] 


VxB—a=0 (1) 


geschrieben worden. Hier ist ® das Tensorfeld der Distorsion (— Deformation + 
Drehung) und @ das Tensorfeld der Versetzungsdichte. Die auch fiir endliche 
Verformungen geltende Formulierung lautet nach Konpo [4] sowie BILBy, 
BuLiLoucH und SmitH [5] (Summationskonvention!) 


Tinie i 4 (On Aj sa 0; Ain) [2 : (2) 


Hier ist J},,4, der in m,/ antisymmetrische Teil der affinen Konnexion J},;;, 
die den natiirlichen Zustand des Kontinuums in den Koordinaten x* des End- 
zustands beschreibt, er reprasentiert die Versetzungsdichte. A,, ist die Trans- 
formation vom natiirlichen in den deformierten Zustand, Aj die dazu reziproke 
Transformation. Die letzte Gleichung entsteht aus der Beziehung 


Le a Ayes a Ay (3) 


durch Antisymmetrisierung beziiglich der Indizes m,/. Wie friiher bemerkt 
({Z], S. 289), ist auch Gl. (3) eine Formulierung des Grundgesetzes. Es ist in 
der Differentialgeometrie wohlbekannt, daB eine Konnexion die Form (3) dann 
und nur dann hat, wenn der zugehérige Riemann-Christoffelsche Kriimmungs- 
tensor [),,,;, verschwindet!. Da wir uns, wie frither begriindet, auf metrische 


' Das Verschwinden des Kriimmungstensors Jj,,,;, in der beschrankten Theorie 
kommt schon bei Bitsy, BuLLouGH und SmrrH vor. 
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Kontinua beschranken, ist J},,,,;, in J, k antisymmetrisch und kann somit voll- 
standig durch den zugehérigen Einstein-Tensor J’ ersetzt werden. Es folgt, 
daB die homogenen Einsteinschen Gleichungen 


(ees! (4) 


eine weitere Formulierung des geometrischen Grundgesetzes in der beschrankten 
Theorie darstellen. 


Hiernach darf geschlossen werden, daB die inhomogenen Einsteinschen Glei- 
chungen 


ese (5) 


in denen B‘’ wegen der Bianchi-Identitat 6 funktionale Freiheitsgrade besitzt, 
eine vollstandige Formulierung des geometrischen Grundgesetzes der allgemeinen 
Theorie (15 Freiheitsgrade) darstellen, wobei die bei der Verformung hervor- 
gerufenen Kontinuitatsst6rungen in der regularen Materie durch die Fremd- 
materie (B‘’) gerade kompensiert werden. 

Wendet man dieselbe SchluBweise auf Gl. (1) an, so kommt man zu einem 
Gesetz 


VxB—a=8, (6) 


in dem der die Fremdmaterie reprasentierende Tensor 6 9 funktionale Freiheits- 
grade zu haben scheint. Ziehen wir von Gl. (6) ihre halbe, mit dem Einheits- 
tensor I multiplizierte Spur (Symbol I) ab, so bleibt unter Beriicksichtigung 
des bekannten Zusammenhangs von Versetzungsdichte und Cosserat-Nyeschem 
Kriimmungstensor K 


Vxp —4,VVx8) 1+ K=5 —451=85. (7) 


Mit der frither bewiesenen Zerlegungsformel ([2], S. 29) 
B=Vs+VxixV+1xO, (8) 
in der s und ) je ein Vektorfeld, t ein symmetrisches Tensorfeld ist, erhalt man 
Vx(VxixV) -OV+K=8'. (9) 


Nun entnehmen wir Gl. (9), daB im Fall 6’=0 und t=0 der Kriimmungstensor K 


die Form a) V hat. Die Realisation erfolgt, wie wir wissen, durch den Angriff 
einer AuBeren Drehmomentendichte. Andere auBere Einfliisse kOnnen nun keines- 


falls eine Kriimmung der gleichen Art, also etwa 9’ V hervorrufen, sie miiBten 
sonst ebenfalls als A4uBere Drehmomente interpretiert werden (wir denken an 
unendlich ausgedehnte Medien). Hieraus folgt, daB K die Form 


K=O0V1®xV7 (10) 


mit einem noch beliebigen Tensorfeld ® hat. Infolgedessen verschwindet die 
linke Seite von Gl. (9) bei rechtsseitiger Divergenzbildung identisch, und es ist 
somit §’- V =0, d.h. der Tensor 6 verliert 3 seiner funktionalen Freiheitsgrade. 
Man iiberzeugt sich leicht, daB andererseits die Gln. (9) bei rechtsseitiger Rotations- 
bildung in die linearisierten Einsteinschen Gleichungen (5) tibergehen (8’ x V= B). 
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Damit ist die Giiltigkeit unserer Behauptung, daB die Einsteinschen Gleichungen 
(5) ein vollstandiger Ausdruck des geometrischen Grundgesetzes sind, auf anderem 
Wege bekraftigt worden. Hiermit werden insbesondere auch die in [/, 6] zuweilen 
gebrauchten Bezeichnungen ,,Grundgleichungen der Eigenspannungsbestimmung 
und ,,Drehmaterie“ iiberfliissig. 

Die Beschrankung von § auf 6 funktionale Freiheitsgrade bedeutet, daB die 
Zahl der Freiheitsgrade aller 4uBeren Einwirkungen nur 12 betragt, im Gegensatz 
zu den statischen und geometrischen ZustandsgréBen, die je 15 Freiheitsgrade 
besitzen. Sind alle a4uBeren Einwirkungen gleich Null, so gibt es infolge der 
besonderen Struktur der Einsteinschen Gleichungen trotzdem von Null ver- 
schiedene Lésungen, es sind dies die zu den Eigenspannungen fiihrenden Lésungen, 
als deren Quellen die Versetzungen auftreten. Die hierzu gehdrigen Freiheits- 
grade sind gerade die 3 iiberzahligen Freiheitsgrade der geometrischen und stati- 
schen GréBen. Offensichtlich handelt es sich hier um Lésungen, die metastabile 
Zustande beschreiben, was in bester Ubereinstimmung mit allem ist, was man 
heute tiber Versetzungen weiB. 

Ich halte es fiir héchst wahrscheinlich, da die Einsteinschen Gleichungen 
iiberhaupt die einzige verniinftige auf den Endzustand bezogene Formulierung 
des geometrischen Grundgesetzes der allgemeinen Theorie durch Differential- 
gleichungen sind, da sie Ableitungen sowohl der Deformationen als auch der 
Kriimmungen enthalten. Dariiber hinaus ist die enge Analogie zur allgemeinen 
Relativitatstheorie sehr ermutigend. Nach den ersten gréBeren Erfolgen zur 
Versetzungsdynamik des Kontinuums von HOLLANDER [7] darf man hoffen, daB 
die Einsteinschen Gleichungen (5) auch bei Hinzunahme der 4. Dimension (der 
Zeit) giiltig bleiben, wobei man jetzt die Indizes in (5) von 1 bis 4 laufen lassen 
sollte. 
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A. Theory of Dilute Suspensions 


G. L. HAND 


Communicated by J. L. ERICKSEN 


1. Introduction 


In this paper the results obtained by JEFFERY [3] for the motion of an 
ellipsoid suspended in a Newtonian fluid are used to formulate a theory of dilute 
suspensions. We show that this theory can be considered a special case of 
ERICKSEN’S theory of anisotropic fluids [/, 2]. 


2. Anisotropic Fluids 

ERICKSEN considers a fluid which is characterized by having a single preferred 
direction at each point, e.g. a suspension of particles of revolution. He lets the 
preferred direction be described by a vector n. ERICKSEN introduces the equations 
o7;—g, where the dot denotes the material derivative. Under the assumptions 
that 

1. the stress tensor ¢;; and g; are functions of the type t,;=1#;;(0, 1;, ;, v; 4) 
8:=8: (0, %, i, U;,;), being linear in m; and v;,, where @ is the density and v, 
is the velocity, 

2.the forms of ¢;; and g; are preserved under all time-dependent proper 
orthogonal coordinate transformations, and 

3. the structure represented by m is symmetric with respect to reflections in 
planes parallel and perpendicular to n, 


he deduces general expressions for ¢;; and g;. Additional assumptions lead to 
simpler forms of ¢;; and g;. These are, that the fluid is incompressible, that 
the particle inertia 9%, is negligible (g;=0), and that m;n;=1. The last assump- 
tion restricts the theory to one representing rigid particles. This gives [2, Eqs. 


(19) and (20) | 
ti, = — 2 O;; + (Ay + Ae dim Me Mn) M4 Nj + 2A dz; + 2A (Ain MN; + 4jn MMi), (1) 
Np, Tad (dp m Nn — dy q Mp Ng a) at Wim Mm > (2) 
where 


(3) 


v5; = 3 (0; 5 a U; é) >» 
5% 5 +25), 4;,=0.- 


The comma denotes differentiation, 6;; is the Kronecker delta, p is an arbitrary 
isotropic pressure, and y and the /’s are constants. The v and A’s might be 
evaluated by comparison with experiment or with some other theory. We in- 
vestigate the latter alternative. 
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3. Jeffery’s Calculation 


JeFFERY considers a rigid ellipsoidal particle suspended in an incompressible 
Newtonian fluid. He assumes that inertial and body forces are negligible and 
that the center of the particle is translated with the mean motion of the fluid. 
JEFFERY assumes that, except in the immediate neighborhood of the particle, 
the flow is steady and homogeneous. Laminar flow prevails throughout the flow 
region, which fills all space exterior to the ellipsoid. 

With the simplified Navier-Stokes and continuity equations, JEFFERY calcu- 
lates the velocity and pressure distributions in the vicinity of the ellipsoid, 
referred to axes coinciding with the principal axes of the ellipsoid. With these 
he is able to calculate the resultant couple exerted on the particle. Since inertial 
forces are assumed negligible, the resultant couple is equated to zero. This 
provides expressions giving the angular velocity w; of the ellipsoid about fixed 
axes which instantaneously coincide with the principal axes of the ellipsoid. 
These expressions are [3, Eq. (37) | 


Ay, O; — Lip Op = Cp R Uk Um, 7» (4) 


where ¢;;,,is the alternating tensor and the a,; are the coefficients of the quadratic 
form describing the ellipsoid. We note that a;;=0 for 1=-7. 

JEFFERY [3, Eq. (56)] calculates the stress on the surface of a sphere centered 
at the suspended particle whose radius R is very large compared to the diameter 
of the particle. The result is approximate since JEFFERY has neglected terms 
of order R-§ and smaller. This stress, referred to axes coinciding with the 


principal axes of the ellipsoid, is 


Lig = Po Oigele ald fer 1061 Rs i 


4%; xj; @ xy ® ; Xj nea 


(5) 


Fee Fe® Re Ré5 
where 
DA (6) 
Gy SiGe dating tras 
6 fo 2 Bo (a? +b?) 2 Bo (a? +b?) 
ihe Gee d,1(Bo —x%) d 
Aes : Sas, 22 | 11\P0 0 ee 
Pq 2 By (a2-+b?) 4b? ag l 12b? BY a 4b? ag (7) 
Beats — azs ds | a, 1 (Bo — 260) 
2 Bo (a? +b?) 4D? xg 4b% a5 ' © 1267 Boas 


ju is the coefficient of viscosity, a, «9, Bo, and fo are given by [3, Eq. (10)], 
and the equation of the ellipsoid of revolution is 

ba de ar 

Sto =. (8) 
We consider an ellipsoid of revolution because the results may be compared 
with formulae occurring in the previously mentioned theory of anisotropic fluids. 
In (5) and (7) the d;; are not the actual rate-of-deformation components of the 
sphere but the constant values approached as R->co. The Ay, remain finite 
as a and b approach zero since each term in A,, is proportional to L%, L being 
some characteristic length of the ellipsoid. t;; becomes the Newtonian stress 


tensor for homogeneous flow of the solvent as the diameter of the suspended 
particle goes to zero. 
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4. Theory of Dilute Suspensions 

We wish to formulate a theory of dilute suspensions based on the results 
of JEFFERY’S calculations. If we restrict the theory to 

1. incompressible fluids, 

2. laminar flow, | 

3. ellipsoidal suspended particles, 

4. a suspension which is very dilute, and 

5. negligible particle inertia, 
then we should be able to apply JEFFERY’s results. 

We must formulate a constitutive equation for the stress. Considering equa- 
tion (5), we have an approximate relation for the stress at a distance from the 
particle which is large compared to the particle diameter. This stress is referred 
to axes coinciding with the principal axes of the ellipsoid. We can average this 
stress over a volume and get the average stress referred to fixed axes which 
instantaneously coincide with these principal axes. We interpret this average 
stress as the macroscopic observable stress of the suspension. We find the space 
average by using the divergence theorem and noting that ¢;; ;=0. 


V S 


where V is the volume over which the average is taken and S its surface*. We 
consider this volume to be a sphere, and, using (5) in (9), we note that the average 
stress becomes exact as the radius of the sphere becomes infinite. Eq. (9) becomes 


Sa ae 4%; 
tp 9 Oe aka Xy+ 10 Apg| oe 0: Fea iP ten hen 
— (% 5 Xp %q Opp + Xi Xp Xp Og5 + %j Xp %q Ops + %; Xp Xp 6,,)|} da; 

1 
= 5 (— 9 6;;+2 d;;) x, da, + 20H Apa. Xz Xp Hq A Ag — oe 
Sana 0°87 a3; ROG; V RS i Xp %q Uap 

S S 
10uA 
S 


since on the boundary of the sphere 9, d;;, Ap, and R are constants. Using 
the divergence theorem, we can integrate and obtain 


44a; = i iy g AV = ~ 0 R yj, 


Ss V 


b*, Ky Xqh Ay, = i (%; Xp Ogn + Xi %q Opn + Xp %q Oin) AV 
S V 


SPORE a 289) (11) 


Re (%; 69444 + %j Opp 44g + Xp O;4 444) 


5 
_ take (Ope ade 0n On 


*x We assume that #,;; ;=0 Ea, the ellipsoid and that ¢;; da; is continuous across 
the surface of the ellipsoid. 


6* 
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Using these results and noting that A,;;=0, we find that (10) has the form 


Py Oe aa a Ay 


8 Loi be 
=— fo Oij + 26 a; + — =F Ag 
where Y is the ratio of particle volume to fluid volume, namely, 
ie sab? ab (13) 


ae REPS 


In JEFFERY’s stress (5) terms of order R-§ and smaller were neglected. In (12) 
the corresponding terms would be those of order W', which are negligible compared 
to terms of order / as long as W is sufficiently small compared to unity. 

Also, we note that we can replace the velocity gradients by their average 
values, since 


1 1 
=> [ %s4 = 7D 42,944, = 04,7 (14) 
V S 


We have averaged as in (9) above, and v,, is considered as being constant in 
regions sufficiently far from the particle. We interpret this average velocity 
gradient as the macroscopic observable velocity gradient which describes the 
motion of the suspension. 

Equation (12) provides the constitutive equation for the stress for a dilute 
suspension of particles which satisfy the restrictions mentioned above, at least 
for cases where the flow of the suspensions is homogeneous. We assume that 
the same equations are applicable also to inhomogeneous flow. 


Now we can show that (12) is just a special case of ERICKSEN’s theory of 
anisotropic fluids. The restrictions presented in §2 give ERICKSEN’s stress the 
simplified form of (1). We choose our coordinate axes in (1) so as to coincide 
with the principal directions of the ellipsoid and interpret ; as a unit vector 
in the direction of the unequal principal axis. Then m,=1, n,=3=0, and (1) 
becomes 


tj = — BOig + (Ay + Ag dys) 0; mj + dg d;; + 224 (dz nj, +5. 4), (15) 


or, In more detail, 


hy =—PtAt (Ag+ 245+ 414) dy, 


beg =— p+ 2A3 doo, 


(16) 
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The mean stresses from (12) are 


iy 4uv 
Np Pe (21 aogat 
2nVv 2nV (By— of 
too =.— fy + (2, aes 7) des 2A ae day, 


7 2uVv t (Be 

ts3 = — pot (2 aie 7) dos + ae a 11> 

e (17) 
4uv 

iy = (24 ab®p ab Ja, 


at a 
A3s= (21 r ab? B , 9 (a2 +b?) } as, 


i 2 us 
isa = (20 + a otag) 420 


For comparison we assume that the d;,’s occurring in (16) and (17) both describe 
the motion of the suspension and that the ¢,; of (16) and ¢;; of (17) both represent 


the stress in the suspension. If we compare to, and es and let As=u(1+0), 
we have 
= (18) 
ey,” 


Comparing ¢,, and ¢,5 or 4,5 and 4,3, we have 


uv 2, 1 
a b2 ( Be (a? +b?) b2 =) : (19) 


Comparing #,, and fy, or Bs and t33, we have 


i Wee 


ist 2uV (xo — Bo) 4 2u¥ (x9 — Bo) 
p Po 3a b4 ag Br = ho + 3a bt ae Cae d;; nN; N;. (20) 
In the theory of anisotropic fluids, # is an arbitrary pressure. Finally comparing 
t,, and ¢,,, we have 


pees ur ( Xx eon ts 4 
2 ab? \ bt ag BY | b2ag = BG (a? +2) /’ 


Ae Oe oe 


Thus either (4) with the constants as just evaluated or (12) may be used as the 
constitutive equation for the stress in this dilute suspension. 

We also wish to know the trajectory of a particle. Since the particle inertia 
is assumed negligible, the center of the particle will be translated with the mean 
motion of the surrounding fluid. The rotation of the particle is given by (2) 
or (4). These give the angular velocities about axes which instantaneously coincide 
with the principal axes of the ellipsoid. To see that (2) and (4) are equivalent, 
we note that the angular velocity occurring in (4) is defined so that %;=e;;,.0; 4. 
We can write w; as the resultant of a vector perpendicular to m,; and one parallel 
to it, the part perpendicular to 1; being given by 


aj} = Cijk 1; Np « (22) 
Then (2) becomes 
Or = Ci7n(V Tap + Wap) Mp Mj. (23) 
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The theory of anisotropic fluids does not account for the component of «, parallel 
to n;. Comparing (23) with (4), we find 
a2— b* 
ane a+b?’ (24) 

where a is the length of the unequal half-axis and } is the length of the equal 
half-axes for the ellipsoid of revolution. Hence y is negative for disk-like particles 
and positive for rod-like particles. Equation (4) or (2) with the constant evaluated 
in (24) provides the trajectory of a suspended ellipsoid of revolution. 

We now have the necessary equations for the laminar motion of a dilute 
suspension of ellipsoids of revolution in an incompressible fluid. The standard 
form of the equations of conservation of mass and momentum give 


d;;=0, 00;=t,;+h, (25) 


where f; is the body force. The stress is found from (12) or (1) with the constants 
evaluated above, and the angular velocity of a particle is provided by (2) or (4). 


This work was supported by Grant NSF-G 11316 from the U.S. National Science 
Foundation. 
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